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Research  during  the  past  year  has  l&d  to  some  significant  results  which  have 
profound  implications  in  the  area  of  vapor  cloud  explosion  in  general,  and  in 
the  area  of  "shockless"  initiation  of  detonation  in  particular.  In  brief  we 
have;  i)  demonstrated  experimentally  that  initiation  of  detonation  can  be 
achieved  by  injecting  a  chemical  catalyst  into  a  fuel-air  mixture,  ii) 
demonstrated  experimentally  that  direct  initiation  can  be  achieved  by  seeding^ 
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an  explosive  mixture  with  free  radicals  of  the  appropriate  concentration  and 
spatial  distribution,  iii)  demonstrated  experimentally  that  direct  initiation 
can  be  achieved  via  turbulent  mixing  between  an  explosive  mixture  and  its 
combustion  products,  iv)  demonstrated  experimentally  that  unconfined  turbulent 
flame  speeds  exceeding  400  m/s  can  readily  be  achieved  in  a  mixture  as  inert 
as  methane  air,  v)  clarified  through  experiments  the  role  of  confinement  on  the 
propagation  of  detonations  near  the  detonability  limits,  and  vi )  assessed  the 
far  field  destructive  potential  of  FAE  weapons  by  calculating  numerically  the 
effective  blast  energy  of  non-ideal  blast  waves  generated  by  vapor  cloud  explo¬ 
sions. 
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Summary  of  Progress 

Our  recent  studies  have  led  to  some  significant  results  which 
have  profound  implications  in  the  area  of  vapor  cloud  explosion  in  gene¬ 
ral,  and  in  the  area  of  "shockless"  initiation  of  detonation  in  particu¬ 
lar.  In  brief  we  have;  i)  demonstrated  experimentally  that  initiation 
of  detonation  can  be  achieved  by  injecting  a  chemical  catalyst  into  a 
fuel-air  mixture,  ii)  demonstrated  experimentally  that  direct  initiation 
can  be  achieved  by  seeding  an  explosive  mixture  with  free  radicals  of  the 
appropriate  concentration  and  spatial  distribution,  iii)  demonstrated  ex¬ 
perimentally  that  direct  initiation  can  be  achieved  via  turbulent  mixing 
between  an  explosive  mixture  and  its  combustion  products,  iv)  demonstra¬ 
ted  experimentally  that  unconfined  turbulent  flame  speeds  exceeding  400  m/s 
can  readily  be  achieved  in  a  mixture  as  inert  as  methane  air,  v) 
clarified  through  experiments  the  role  of  confinement  on  the  propagation 
of  detonations  near  the  detonability  limits,  and  vi)  assessed  the  far  field 
destructive  potential  of  FAE  weapons  by  calculating  numerically  the  effec¬ 
tive  blast  energy  of  non-ideal  blast  waves  generated  by  vapor  cloud  explo¬ 
sions. 


During  the  past  year  a  significant  milestone  towards  the  realiza¬ 
tion  of  an  FAE  III  device  has  been  reached  at  McGill  in  that  detonation  in 
a  premixed  fuel-air  system  (propane-air)  solely  by  injecting  a  chemical 
agent  (fluorine)  has  been  achieved  on  a  laboratory  scale.  Laboratory  scale 
experiments  are  of  necessity  small  scale  and  confined,  and  a  considerable 
research  effort  has  been  made  to  identify  the  fundamental  underlying  mecha¬ 
nisms  and  to  establish  quantitative  criteria  for  the  scaling  up  of  these 
experiments  in  preparation  for  large  scale  bag  tests  to  be  performed  at 
Eg! in  Air  Force  Base.  A  review  of  the  progress  in  the  area  of  Pyrophoric 
Initiation  of  Detonation  in  Fuel-Air  Mixtures  including  a  discussion  of 
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criteria  for  scale-up  and  removal  of  confinement,  as  well  as  a  discussion 
of  the  dispersion  and  mixing  requirements  are  given  in  Appendix  I. 

The  role  of  both  large  and  small  scale  turbulence  in  augmenting 
the  explosion  of  a  fuel-air  cloud  have  been  further  clarified  by  our  inves¬ 
tigations  of  turbulent  flame  propagation  and  acceleration  in  the  presence 
of  obstacles.  These  investigations  have  established  that  extremely  high 
flame  speeds  and  large  overpressures  can  be  obtained  through  the  generation 
of  large  scale  turbulent  eddies  and  small  scale  turbulence  in  the  unburned 
gas  by  placing  repeated  obstacles  in  the  flame  path.  Laboratory  scale  ex¬ 
periments  on  unconfined  cylindrically  expanding  methane-air  flames  with 
repeated  obstacles  in  the  flame  path  using  both  streak  and  schlieren  diag¬ 
nostic  observations,  together  with  larger  scale  experiments  performed  in 
collaboration  with  Professor  Wagner  at  the  University  of  Gdttingen,  have 
deMnonstrated  that  flame  speeds  in  excess  of  400  m/s  with  associated  over¬ 
pressures  of  the  order  of  0.65  bar  can  be  generated  in  an  obstacle  environ¬ 
ment.  A  simple  theoretical  model  of  the  flame  acceleration  process  based 
on  a  positive  feedback  mechanism  between  the  flow  field  generated  by  the 
flame  propagation  over  obstacles,  and  the  burn-out  of  the  large  turbulent 
eddies  predicts  an  exponential  increase  in  flame  velocity  which  is  in  good 
agreement  with  the  experimental  results.  The  current  status  of  our  inves¬ 
tigations  on  turbulent  flame  propagation  including  detailed  discussions  of 
experimental  and  theoretical  results  is  given  in  Appendix  II. 

The' detonabil i ty  limits  determine  the  conditions  beyond  which  an 
explosive  mixture  can  no  longer  support  a  detonation.  At  the. present  time 
it  is  not  possible  to  predict  the  detonabil ity  limits  of  a  given  mixture 
theoretically  and  the  limits  must  be  determined  experimentally.  Experiments 
to  determine  these  limits  are  usually  performed  in  detonation  tubes,  how¬ 
ever  the  influence  of  confinement  and  the  initiation  source  on  the  propagation 


of  detonations  is  not  understood.  Our  investigation  of  the  propagation  of 
detonations  in  long  tubes  (14m)  of  diameter  4.8  cm  and  14.5  cm  near  the 
lean  limit  of  detonability  for  ethylene-air  at  1  atm  shows  that  the  propa¬ 
gation  of  marginal  detonations  are  strongly  influenced  by  both  the  initial 
conditions  and  the  boundary  conditions.  In  fact,  completely  different 
phenomena  are  observed  for  the  same  mixture  composition  depending  on  the 
initial  and  boundary  conditions.  Transverse  spinning  waves  with  wavelengths 
characteristic  of  the  tube  rather  than  the  mixture  are  found  to  play  a  domi¬ 
nant  role  in  the  propagation  of  near-limit  detonation  waves.  Because  of 
this  strong  influence  of  the  confinement  a  criteria  for  characterizing  the 
detonability  limits  for  unconfined  detonations  based  on  the  onset  of  a 
single  head  spin  has  been  proposed.  Associated  with  this  limit  there  is 
then  a  characteristic  length  equal  to  the  tube  diameter  (or  spin  pitch), 
and  by  extrapolating  from  results  in  various  diameter  tubes  the  limit  for 
unconfined  detonations  in  a  cloud  of  given  size  can  be  determined  once  the 
minimum  number  of  transverse  waves  required  for  unconfined  detonation  is 
determined.  A  detailed  discussion  of  the  Influence  of  Confinement  on  the 
Propagation  of  Detonations  near  the  Detonability  Limits  based  on  our  results 
is  given  in  Appendix  III. 

The  far  field  destructive  potential  of  fuel-air  explosive  weapons 
has  been  assessed  by  calculating  numerically  the  effective  blast  energy  of 
the  non-ideal  blast  waves  generate  by  vapor  cloud  explosives.  For  hydro¬ 
carbon-air  and  oxygen  explosives,  the  effective  blast  energies  represent 
approximately  37  and  2 Tl  of  the  respective  combustion  energy  of  the  mixture, 
irrespective  of  the  mode  of  combustion  (i.e.,  detonation  or  volumetric  ex¬ 
plosion).  Based  on  the  effective  blast  energy,  far  field  equivalence  of 
non-ideal  blast  waves  from  fuel -air  and  oxygen  explosives  is  demonstrated. 

In  other  words,  the  far  field,  non-ideal  blast  parameters  (i.e.,  overpressure. 
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static  and  dynamic  impulses)  scale  according  to  the  explosion  length  based 
on  the  effective  blast  energy  and  can  be  radily  calculated  from  a  single 
chart  to  assess  the  destructive  potential  of  fuel-air  and  oxygen  explosive 
weapons.  Although  all  ideal  point  blast  waves  from  solid  explosive  weapons 
are  equivalent  and  the  corresponding  blast  parameters  scale  according  to 
the  explosion  length  based  on  the  total  combustion  energy  of  the  charge,  the 
present  study  indicates  that  ideal  and  non-ideal  blast  waves  are  not  equiva¬ 
lent  in  the  far  field.  In  other  words,  at  the  same  scaled  distance,  far 
field  blast  parameters  from  ideal  and  non-ideal  blast  waves  are  no  equal. 
Based  on  the  matching  of  the  overpressure  curves  in  the  far  field,  the  equi¬ 
valent  point  source  energy  for  the  non-ideal  blast  wave  from  a  fuel -air 
and  oxygen  explosive  weapon  is  three  times  greater  than  the  energy  for  a 
point  blast.  Therefore,  the  far  field  destructive  potential  of  fuel -air 
and  oxygen  explosive  weapons  is  found  to  be  superior  to  that  of  solid  ex¬ 
plosive  weapons.  A  detailed  discussion  and  presentation  of  the  key  results 
of  this  investigation  on  blast  waves  from  fuel-air  explosives  are  given  in 
Appendix  IV. 

In  summary  our  investigations  on  the  transition  to  detonation, 
the  influence  of  confinement  on  the  detonability  limits  and  flame  accelera¬ 
tion  have  provided  and  are  continuing  to  provide  important  criteria  and 
guidelines  for  both  accidental  explosions  and  the  development  of  an  FAE  III 
device.  Our  understanding  of  these  phenomena  are  now  at  the  stage  where 
critical  calculations  and  experiments  are  being  performed.  In  the  case  of 
FAE  III  for  example,  large  scale  feasibility  experiments  are  being  plannejl— 
this  year. 


1 .  Pyrophoric  Initiation  of  Detonation  in  Fuel-Air  Mixtures 

The  present  is  a  summary  report  of  the  progress  achieved  during 
the  fiscal  year  of  1979  in  the  investigation  of  the  mechanisms  of  initiation 
of  detonation  in  fuel-air  mixtures  via  a  chemical  catalyst.  This  work  has 
direct  relevance  to  the  development  of  the  single  event  FAC  III  device.  The 
main  objective  of  our  work  at  McGill  University  is  a  continuing  process  of 
research  to  demonstrate  the  feasibility  of  this  mode  of  initiation  for  the 
actual  realization  of  an  FAE  III  device  based  on  the  co-dispersal  of  a  fuel 
and  a  chemical  catalyst  into  air.  The  approach  in  the  present  program  is 
to  identify  the  fundamental  underlying  mechanisms  and  to  establish  quanti¬ 
tative  criteria  for  this  mode  of  initiation  on  a  laboratory  scale  to  be  fol¬ 
lowed  up  by  feasibility  demonstration  studies  in  large  scale  bag  tests  in 
collaboration  with  and  at  Eglin  Air  Force  Base. 

In  early  1979  experiments  using  fluorine  injection  into  heavy 
hydrocarbon-air  mixtures  (butane-air,  propane-air)  were  commenced.  By  nrid- 
1979,  direct  initiation  of  detonation  in  a  premixed  fuel-air  system  (pro¬ 
pane-air)  solely  by  injecting  a  chemical  catalyst  (fluorine)  was  achieved. 

A  schematic  diagram  of  the  apparatus  is  shown  in  Fig.  1.  In  this  case 
charges  of  fluorine  and  oxygen  (optimum  ratio  1 07.'  fluorine  and  90«  oxygen) 
were  counter-injected  against  a  small  charge  of  propane  into  a  premixed 
reservoir  of  stoichiometric  propane-air  mixture.  The  fluorine  interacts 
with  the  propane  to  form  an  unstable  radical  (CJU).  The  oxidation  reactions 
then  progress. rapidly  in  the  presence  of  oxygen  via  a  chain-breaking  mecha¬ 
nism  leading  to  ignition  in  the  propane-air  mixture  and  ultimately  initia¬ 
tion  of  detoation.  No  detonation  initiation  was  achieved  in  other  hydro¬ 
carbon-air  mixtures  except  propane-air  even  when  these  were  more  detonable 
than  propane-air  (eg.,  ethylene-air).  This  confirms  the  effectiveness  of 
the  fluorine-heavy  hydrocarbon  interaction  in  radical  production  as  indicated 
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by  von  Elbe  originally.  The  diagnostics  in  the  present  experiments  include 
a  photodiode  to  detect  the  onset  of  chemiluminescent  reactions*  the  signal 
from  which  is  used  to  trigger  the  oscilloscope.  Two  pressure  transducers 
are  used  to  monitor  the  onset  of  detonation:  one  in  the  vicinity  of  injec¬ 
tion  either  above  or  to  the  side,  the  other  at  the  end  of  the  extension  tube 
as  illustrated  in  the  schematic.  A  typical  pressure  record  illustrating 
the  case  where  detonation  level  pressure  was  achieved  at  the  end  transducer 
(upper  beam)  is  displayed  in  the  oscilloscope  trace  in  Fig.  2.  The  lower 
beam  illustrates  the  pressure  record  in  the  vicinity  of  the  injection.  To 
be  noted  is  the  fact  that  it  arises  later  than  the  pressure  spike  at  the 
extension  tube  end  indicating  that  the  development  of  detonation  was  in 
the  region  between  the  injection  and  the  end  of  the  extension  tube.  This 
is  compatible  with  the  SWACER  mechanism  concept  of  detonation  iriitation. 

Thus  the  feasibility  of  catalytic  detonation  initiation  in  a  fuel-air  mix¬ 
ture  has  been  demonstrated  on  the  laboratory  scale. 


2.  Criteria  for  Deconfinement  and  Scale-up 


Laboratory  scale  experiments  of  practical  necessity  must  be  small 
scale  and  confined.  To  demonstrate  the  feasibility  of  an  FAE  III  single 
event  device  conclusively,  large  scale  unconfined  experiments  must  be  ulti¬ 
mately  performed.  In  the  present  program  it  is  envisioned  that  these  ex¬ 
periments  will  be  carried  out  at  Eglin  Air  Force  Base  using  bag  tests.  To 
progress  from  laboratory  scale  to  bag  test  scale,  effects  of  deconfinement 
and  scale-up  must  be  considered  parametrically.  It  is  well  known  that  deto 
nations  in  confined  tubes  can  be  readily  achieved  even  in  the  most  insen¬ 
sitive  fuel-air  mixtures  (eg.,  methane-air)  as  demonstrated  by  Gerstein  and 
Kogarko.  However,  detonation  initiation  under  unconfined  conditions  is 
ve>'y  difficult  in  methane-air,  requiring  large  initiating  charges  of  solid 
explosive  (viz.  Bull  and  Benedick).  Thus  confinement  helps  detonation 
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initiation  and  propagation.  For  example,  it  is  well  known  that  if  a  detona¬ 
tion  is  initiated  and  propagates  as  a  well  established  wave  in  a  confined 
tube,  it  may  fail  when  it  emerges  abruptly  into  an  unconfined  geometry  unless 
the  diameter  of  the  tube  from  which  the  detonation  emerged  has  a  diameter 
greater  than  some  critical  tube  diameter  for  the  particular  mixture.  This 
critical  tube  diameter  for  unimpeded  transmission  of  a  detonation  from  con¬ 
fined  linear  propagation  to  a  fully  unconfined  geometry  has  been  found  to 
be  of  great  fundamental  significance.  It  can  be  related  to  the  critical 
energy  for  direct  initiation  of  detonation  (viz.  Lee  and  Matsui)  and  is  a 
measure  of  the  geometrical  detonability  limit  which  indicates  the  minimum 
cross-sectional  dimension  or  thickness  of  a  fuel -air  cloud  below  which  a  de¬ 
tonation  can  no  longer  sustain  its  steady  propagation.  A  table  of  experi¬ 
mentally  established  measurements  and  estimates  of  the  critical  tube  diameter 
for  a  range  of  fuel -oxidizer  mixtures  ranging  from  fuel-oxygon  to  fuel-air 
are  presented  in  Fig.  3.  These  results  indicate  that  the  critical  tube  dia¬ 
meter  for  propane-air  would  be  in  excess  of  two  meters.  This  would  infer 
that  no  detonation  propagation  would  be  possible  in  a  propane-air  cloud  under 
fully  unconfined  conditions  if  the  cloud  thickness  were  less  than  this. 
However,  in  the  case  of  a  bag  test  this  critical  dimension  would  be  somewhat 
less  because  of  the  partial  confinement  due  to  ground  effect.  In  the  case 
of  a  large  scale  FAE  feasibility  experiment  the  implication  is  that  one  would 
have  to  generate  a  region  of  pyrophorically  initiated  chemical  activity  which 
would  result  in  an  initiating  detonation  wave  extending  over  a  dimension  at 
least  of  the  order  of  the  critical  tube  diameter  for  that  particular  mixture. 
Moreover,  the  energetics  of  the  initiation  process  would  further  require 
that  the  energy  contained  in  the  initiating  volume  should  be  at  least  of 
the  order  of  the  critical  energy  for  direct  initiation  of  detonation  in  the 
mixture.  These  considerations,  therefore,  provide  quantitative  guidelines 
whereby  the  critical  volume  of  fuel -air  mixture  which  must  be  pyrophorically 
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activated  for  detonation, initiation  can  be  estimated.  It  is  clear  that  the 
initiating  volume  for  propane-air  appears  to  be  prohibitively  large  (i.e., 
a  cross-sectional  dimension  in  excess  of  2  meters).  This  would  imply  that 
the  catalyst  (fluorine-oxygen)  would  have  to  be  dispersed  simultaneously 
over  a  cross-section  in  excess  of  2  meters  and  perhaps  a  half  meter  deep 
in  the  propane-air  cloud.  Practical  problems  of  the  injection  process  would 
therefore  appear  to  preclude  propane-air  as  the  FAE  cloud.  More  reasonable 
dimensions  v;ould  indicate  that  ethylene-air  can  sustain  detonation  in  a 
dimension  of  cloud  less  than  one  me^ter  (i.e.,  based  on  a  critical  tube  dia¬ 
meter  of  80  cm).  Although  no  critical  tube  diameter  data  exist  for  MAPP- 
air  or  propylene  oxide-air  mixtures,  critical  tube  diameter  values  for  these 
mixtures  can  be  estimated  via  a  simple  analysis  (viz.  Lee  and  Matsui)  from 
the  critical  initiation  energy  values,  say  those  obtained  by  Parsons  et  al . 
Based  on  such  estimates  MAPP-air  would  have  a  critical  tube  diameter  approxi¬ 
mately  that  of  ethylene-air  (i.e.,  80  cm)  while  propylene  oxide-air  would 
be  even  less  than  that.  In  terms  of  an  FAE  bag  test  we  are  therefore  talk¬ 
ing  of  a  bag  of  cross-section  slightly  in  excess  cf  a  meter  and  several 
meters  long  with  a  pyrophoric  injection  cross-section  of  the  order  of  50-80  cm. 
It  is  important  to  recall  that  since  the  active  ingredients  for  catalytic 
initiation  are  propane,  fluorine  and  oxygen,  these  constituents  must  be 
provided  in  the  initiating  region  of  the  fuel-air  vapor  cloud  which  may 
now  be  either  ethylene-air,  MAPP-air  or  propylene  oxide-a:r.  It  is  obvious 
that  for  the  initiation  process  to  be  successful  the  dispersion  time  must 
be  shorter  than  the  catalytic  chemical  activation  time  which  is  typically 
short,  being  of  the  order  of  a  few  milliseconds  in  this  case.  .  Consequently, 
in  the  latter  part  of  1979,  the  main  efforts  in  the  FAE  program  at  McGill 
have  concentrated  on  the  study  of  rapid  simultaneous  injection  and  dispersal 
techniques  via  multipoint  sources  spread  over  a  large  cross-sectional  area 
in  anticipation  of  the  approach  to  be  used  in  the  bag  test  at  Eglin  Air  Force 
Base. 
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3.  Rapid  Injectjon, sponsion  and _  Mjxjncj_o_f_  Gasps 

In  the  original  small  scale  catalytic  initiation  experiments  at 
McGill,  manually  activated  toggle  valves  wore  used  successfully  for  the 
injection  of  the  pyrophoric  constituents  into  the  fuel-air  mixture.  This 
proved  adequate  because  only  very  small  regions  of  chemical  activation 
wore  involved  and  the  dispersion  was  sufficiently  rapid.  In  the  uncon¬ 
fined  large  scale  bag  test  context,  multipoint  simultaneous  injection  and 
dispersal  on  a  millisecond  time  scale  would  be  necessary  for  reasons  al¬ 
ready  outlined.  Clearly  the  use  of  remotely  activated  mechanical  toggle 
valves  with  their  relatively  slow  opening  Limes  and  synchronization  prob¬ 
lems  would  preclude  their  use  in  this  case.  The  requirement  for  detonation 
initiation  in  a  bag  test  would  be  for  very  rapid  injection,  deep  penetra¬ 
tion,  broad  lateral  spread  and  intense  turbulent,  mixing  of  the  catalytic 
ingredients  with  the  fuel-air  cloud  simultaneously  from  several  injection 
ports  to  cover  the  critical  cross-sectional  dimension  equivalent  to  the 
critical  tube  diameter  all  on  a  millisecond  tinu*  scale.  Such  requirements 
would  indicate  an  active  rupture  type  dispersion  system  and  experiments 
were  initiated  to  study  injection  of  gas  charges  into  an  ambient  medium 
via  a  piston  driver-diaphragm  scheme  (Fig.  4).  Argon  gas  which  has  vir¬ 
tually  the  same  atomic  weight  as  the  fluorine  molecular  weight  was  used 
initially  as  the  injected  gas  into  ambient  air  lor  the  penetration  and 
dispersion  studies.  A  typical  spark  seblieren  sequence  of  photographs 
illustrates  a  typical  injection  sequence  in  I  ig.  b,  Sufficiently  rapid 
injection  and  adequate  penetration  was  achieved  in  this  way.  I  ig.  6  dis¬ 
plays  typical  transient  shock-jet  trajectories  to  confirm  this.  However, 
when  fluorine  was  subsequently  used,  its  corrosive  properties  affected  the 
piston  sliding  trajectory  I  rom  shot  to  shot  in  view  of  Lh<>  fight,  clearances 
in  the  honed  cylinders  required  for  sealing  purpose.,  consequently  this 
scheme  of  injection  was  deemed  too  intricate  and  unreliable  for  field  test 
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situations.  An  alternate  scheme  of  simultaneous  impulsive  rupturing  of 
an  array  of  prepressurized  glass  tubes  was  considered  more  promising. 

I. I.  Glass  of  UTIAS  has  done  considerable  work  on  the  gasdynamic  flow 
field  generated  from  the  rupture  of  cylindrical  and  spherical  pressurized 
glass  vessels  in  the  early  sixties.  Of  particular  interest  here  is  the 
very  turbulent  spread  of  the  contact  surface  between  the  pressurized  and 
ambient  gases  as  a  result  of  the  obstructing  effect  of  the  shattered 
fragments  of  the  glass  containers  which  are  initially  immobilized  because 
of  their  large  inertia.  The  sequence  of  schlieren  photographs  in  Fig.  7 
taken  from  Glass's  book  illustrates  this.  Calculations  at  McGill  using 
a  gasdynamic  numerical  code  have  been  carried  out  for  spherical  charges 
and  compared  with  experimental  streak  photographs  of  Glass  for  identical 
initial  conditions.  The  agreement  between  theory  and  experiment  as  indi¬ 
cated  in  Fig.  8  for  the  motion  of  the  contact  surface  is  remarkable.  The 
particular  example  calculated  here  is  for  a  practically  realistic  situa¬ 
tion  in  the  context  of  an  actual  bag  test.  We  are  contemplating  glass 
flask  charges  of  pyrophoric  gases  with  gas  overpressure  ratios  of  22-25 
(i.e.,  325-350  psig).  As  indicated  in  Fig.  8,  the  maximum  contact  surface 
displacement  for  such  overpressure  would  extend  in  excess  of  three  times 
the  initial  dimension.  For  an  initial  spherical  flask  dimension  of  20-25 
cm  diameter,  the  gas  would  disperse  in  a  time  of  the  order  of  500  wsec 
over  a  spherica1  volume  in  excess  of  75  cm  diameter.  An  added  advantage 
of  this  mode  of  dispersion  apart  from  its  speed  is  the  fact  that  the  iner¬ 
tia  of  the  shattered  glass  fragments  would  provide  a  "screen"  effect 
through  which  the  contact  front  would  penetrate  and  thus  become  highly 
turbulent.  This  would  greatly  enhance  the  rapid  mixing  capability  of  the 
pyrophoric  materials.  In  an  actual  bag  test  initiation  experiment  it  is 
contemplated  to  use  three  flasks  spaced  about  two  flask  diameters  apart. 
The  central  flask  would  be  pressurized  with  propane;  the  other  two  glass 
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sphere  flasks  would  contain  a  pressurized  fluorine-oxygen  mixture  in  the 
ratio  of  1:9  fluorine  to  oxygen.  Fig.  9  indicates  the  sequence  of  events 
in  the  initiation  scheme.  The  central  flask  containing  propane  would  be 
ruptured  first  via  mechanical  fracture  using  mild  detonating  fuse  (MDF). 
The  two  peripheral  flasks  containing  the  fluorine-oxygen  charges  would 
then  be  ruptured  in  the  same  way  but  about  500  psec  later  to  permit  the 
propane  from  the  first  flask  to  have  dispersed  to  its  maximum  extent.  The 
subsequent  rapid  turbulent  mixing  of  the  pyrophoric  materials  over  a  cross 
sectional  dimension  equivalent  to  the  critical  tube  diameter  for  ethylene- 
air,  MAPP-air  or  propylene-oxide-air  mixture  in  the  bag  would  lead  to  deto 
nation  initiation. 
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APPENDIX  II 

Turbulent  Flame  Propagation  and  Acceleration 
in  the  Presence  of  Obstacles* 


*  Paper  submitted  to  the  Proceedings  of  the  Vllth  International  Colloquium 
on  Gasdynamics  of  Explosions  and  Reactive  Systems. 
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**1 .  IntrodiKitioji 

V 

Violent  flame  acceleration  ami  enhanced  transition  to  detonation 

when  turbulence  producing  obstacles  arc  placed  in  the  path  of  a  flame  have 

1  -3 

been  observed  by  various  authors.  However,  at  the  present  time,  the  ques¬ 
tions  of  maximum  flame  speed  and  flame  acceleration  in  unconfined  or  par¬ 
tially  confined  situations  with  turbulence  producing  obstacles  and/or  initial 
turbulence  are  still  far  from  being  resolved,  from  the  practical  point  of 
view  it  is  important  to  determine  the  flame  acceleration  and  maximum  flame 
speed  for  a  wide  range  of  turbulence  conditions  since  the  turbulence  in  a 
cloud  of  combustible  mixture,  which  may  result  from  an  accidental  spill  of 
explosive  gases  and  liquids,  can  vary  widely  depending  on  the  nature  and 
location  of  the  spill,  Cven  if  the  turbulence 'level  is  lev/  prior  to  ignition, 
the  presence  of  obstacles  in  the  form  of  buildings,  vehicles,  pipes,  etc.,  can 
lead  to  large  flame  speeds  and  potentially  dangerous  blast  waves. 

At  the  present  time,  flame  laws  relating  the  flame  speed  of  a  freely 
propagating  flame  to  the  turbulence  field  or  obs  tao  1  o  <;onf  i  gura  tion  arc  not 
available.  Such  flame  laws  are  required  as  input,  to  existing  numerical  gas- 
dynamic  calculations*1  ^  to  evaluate  the  resulting  pressure  waves,  so  that 
the  blast  effects  due  to  a  fuel-air  explosion  in  a  turbulent  environment  can 
be  predicted.  Wo  have  undertaken  a  program  whose  aim  is  t.n  develop  such 
flame  laws. 

In  a  previous  paper*7  we  reported  on  the  dramatic  acceleration  ol 
cyl Indrical ly  expanding  methane-air  flames  when  repented  obstacles  were  placed 
in  the  flame  path.  With  repeated  obstacles  of  the  appropriate  sizes  and 
separations,  flame  speeds  up  to  130  m/s  were  obtained.  Ibis  was  approximately 
24  times  the  speed  observed  with  no  obstacles.  T  he  mechanism  proposed  for  the 

rapid  acceleration  of  the  flame  from  an  initial  flame  speed  of  about  3.3  m/s 

/ 

to  speeds  in  excess  of  100  m/s  over  a  distance  oi  3r..‘»  cm  was  based  on  the 


‘positive  feedback  coupling  between  the  flame  and  the  upstream  flow  produced 
by  the  specific  volume  increase  across  the  flame.  With  repeated  obstacles 
in  the  flame  path,  gradients  in  the  upstream  flow  field  produce  "flame  fold¬ 
ing"  which  increases  flame  area  and  hence  the  rate  of  burning,  leading  to  a 
larger  flow  velocity  and  stronger  flow  field  gradients,  which  then  further 
increases  the  rate  of  burning  and  so  on.  This  positive  feedback  mechanism 
can  thus  lead  to  a  flame  which,  without  some  damping  or  limiting  factors, 
will  continue  to  accelerate  as  long  as  the  obstacle  configuration  remains 
the  same. 

Most  of  the  results  reported  in  deference  1  were  based  on  experi¬ 
ments  in  a  chamber  which  provided  only  30.5  cm  of  flame  travel.  A  few  ex¬ 
periments  wore  performed  in  a  larger  chamber  providing  hi  cm  of  flame  travel. 
These  latter  experiments  confirmed  the  dramatic  influence  of  obstacles  on 
the  flame  propagation,  but  the  results  were  inconclusive  on  the  questions  of 
continued  flame  acceleration  and  maximum  flame  speed.  To  investigate  these 
questions  further  we  constructed  a  larger  chamber  so  that,  cylindrical  flame 
propagation  could  be  observed  over  a  distance  of  120  cm.  Ibis  paper  reports 
on  the  results  of  this  investigation.  Based  on  tnese  results  together  with 
spark  schlieren  photographs  of  the  flame  propagation  ov<w  obstacles  in  a  rec¬ 
tangular  channel,  a  simple  model  of  the  flame  acceleration  process  is  pro¬ 
posed.  This  model  predicts  an  exponential  increase  in  flame  speed  with  dis¬ 
tance  propagated  over  obstacles  which  is  in  good  agreement  with  the  experi¬ 
mental  results. 

2.  experimental  Details 

•  The  investigation  of  flame  propagation  over  obstacles  was  performed 
in  a  2.5m  x  2.5m  chamber  consisting  of  two  wooden  parallel  plates  separated 
by  spacers  at  the  corners.  The  plate  separation  could  be  adjusted  by  changing 
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‘the  spacers.  The  chamber  was  scaled  at  the  periphery  by  plastic  hands  while 
it  was  filled  wit.lt  the  appropriate  gas  mixture.  A  spring  mechanism  was  used 
to  remove  these  plastic  bands  prior  to  ignition  so  that  the  flame  could  ex¬ 
pand  freely  in  the  radial  direction.  The  gas  mixture  of  stoichiometric  me¬ 
thane-air  was  prepared  by  continuous  flow,  and  the  chamber  was  purged  at  least 
10  times  in  order  to  ensure  a  good  mixture  in  the  chamber.  Ignition  was 
achieved  by  a  fast  burning  wire  at  the  center  of  the  chamber,  and  the  subse¬ 
quent  radial  propagation  of  the  flame  was  recorded  by  a  rotating  drum  streak 
camera.  A  schematic  diagram  of  the  experimental  set-up  is  shown  in  figure  1. 

The  obstacle  configurations  wore  identical  to  those  used  in  the 
laboratory  experiments  reported  in  Reference  7,  with  obstacles  consisting  of 
spirals  of  copper  or  plastic  tubing.  The  spiral  obstacles  were  placed  on 
the  top  or  bottom  plate  so  that  a  radially  expanding  flame  would  encounter  a 
series  of  obstacles,  corresponding  to  the  windings  of  the  spirals,  which  ex¬ 
tended  all  the  way  to  the  periphery  of  the  chamber.  A  cross-sectional  slice 
through  the  center  of  the  chamber  illustrating  the  obstacle  configurations  is 
shown  in  Figure  2.  The  obstacle  height  II,  the  pitch  f  and  the  plate  separa¬ 
tion  D  are  also  defined  in  this  figure,  for  these  particular  experiments 
two  types  of  obstacle  spirals  were  used;  i)  a  spiral  of  1.25  cm  diameter  cop¬ 
per  tubing  wi th  a  pitch  of  3.75  cm,  and  ii)  a  spiral  of  A  cm  plastic  tubing 
with  a  pitch  of  10  cm.  The  plate  separation  was  varied  bi-tween  3.5  cm  and 
20  cm. 

To  investigate  in  more  detail  the  mechanism  of  flame  propagation 
over  obstacles  wo  have  also  taken  schl ieren  photographs  of  a  methane-air 
flame  propagating  over  obstacles  in  a  rectangular  channel  of  length  30  cm. 
This  channel  was  constructed  by  placing  wooden  bars  in  the  30.5  cm  radius 
chamber  described  in  Reference  7,  and  the  schl ieren  photographs  were  taken 
through  windows  in  the  two  plates  of  the  chamber.  A  schematic  diagram  of 


'the  channel  with  obstacles  is  shown  in  Figure  3,  where  the  relevant  parameters ; 
pitch  (P),  obstacle  diameter  (ll*)  and  channel  width  (W)  are  also  defined.  The 
height  of  the  channel  was  2 .  iid  cm.  for  some  of  the  experiments  in  the  channel 
one  or  more  of  the  obstacles  were  filled  with  dry  ice.  Three  small  holes  <n 
these  obstacles  allowed  CO^  to  escape.  The  escaping  CO^  provides  an  excellent 

schlieren  tracer  for  the  flow  field  ahead  of  the  flame. 

* 

3 •  General  Cjon_s_idej'_a_ti_oji_s 

As  discussed  in  Reference  7,  the  dramatic  influence  of  repeated  ob¬ 
stacles  on  the  flame  propagation  can  be  understood  in  terms  of  the  positive 
feedback  coupling  between  the  flame  itself  and  the  flow  ahead  of  the  flame. 

When  obstacles  are  placed  in  the  flame  path,  the  upstream  flow  field  will  be 
that  characteristic  of  the  obstacle  configuration.  for  repeated  obstacles 
placed  on  the  bottom  and/or  top  plates  of  the  chamber,  the  main  features  of 
the  flow  ahead  of  the  flame  are  a  standing  eddy  in  the  wake  region  behind  the 
obstacles  which  is  separated  from  the  outer  flow  region  by  a  shear  layer.  T ho 
flow  pattern  ahead  of  the  flame  and  the  entrainment  into  the  standing  eddy 
for  low  flame  velocities  is  illustrated  by  the  spark- sc h 1 i eron  photograph 
shown  in  Figure  4a.  In  this  case  the  second  and  the  fourth  obstacles  have 
been  filled  with  dry  ice  and  the  resulting  CO.,  tracer  clearly  shows  a  stand¬ 
ing  eddy  region  between  the  second  and  third  obstacle.  As  the  flame  encounters 
this  flow  field  it  will  become  "folded"  due  to  the  gradient,  in  the  flow  field, 
leading  to  a  flame  consisting  of  a  curved  leading  flame  front  in  the  outer 
flow  region  with  a  trailing  flame  in  the  standing  eddy  between  the  obstacles. 
This  Is  illustrated  in  the  spark  schlieren  photograph  shown  in  Figure  4h.  This 
photograph  also  shows  the  entrainment  of  the  flame  eono  into  the  eddy.  Tor 
larger  flame  speeds  and  thus  larger  flow  velocities  ahead  of  the  flame,  the 
shear  layers  and  the  wakes  of  the  obstacles  will  become  turbulent  and  as  the 
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velocity  increases  the  turbulent  flow  region  will  encompass  more  a  ml  more  of 
the  flow  field.  This  is  illustrated  in  figure  ‘>u  whore  dry  ice  in  the  fourth 
obstacle  provides  the  COp  tracer.  The  turbulent,  wrinkling  and  entrainment  of 
the  flame  zone  into  the  wake  region  can  also  be  seen  in  this  figure.  The 
smaller  scale  turbulent  wrinkling  of  the  flame  is  further  illustrated  by  the 
spark-schl  i cron  photograph  shown  in  Figure  f>b.  Tltese  spark-schl  i eron  photo¬ 
graphs  confirm  the  picture  of  flame  propagation  over  obstacles  proposed  in 

Reference  7.  A  schematic  diagram  illustrating  Litis  physical  picture  is  shown 

♦ 

in  Figure  6. 

As  discussed  in  Reforone  7,  the  turbulent,  burning  velocity  Sy  based 
on  the  average  heat  release  rate  per  unit  projected  or  frontal  f'-  me  area  is 
given  by: 

Sj'i'.roS*-'-  " 


whore  the  integral  is  over  the  area  of  the  "flam*'  fold"  shown  in  Figure  A,  f> 
is  the  local  burning  velocity  at  the  "flame  fold"  surface,  11  is  the  plate 
separation  and  r  is  the  radial  position  of  the  I  lame  trout.  Assuming  zero 
flow  velocity  in  the  onhorned  gases  behind  the  flame,  the  displacement  velocity 
of  t ho  unburned  gases  just  ahead  of  the  flame  averaged  over  the  plate  separa¬ 
tion  i.s  given  by: 


<v-  •  »s, 


(2) 


where  o  and  are  the  average  densities  of  the  unburned  gas  ahead  of  t ho 


flame  and  the  burned  gas  behind  the  flame  zone,  respectively.  Thus,  assuming 
that  the  advance  of  ttie  flame  is  uniform  across  the  plate  separation  and  that 
transient  effects  can  be  neglected,  the  corresponding  flame  speed  is  given  by 


Rf  =  ~u  ST 
f  **b  T 


(3 
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'  lor  low  f  1  .imo  speeds  in  the  initial  ol  the  flame  acceleration, 

the  flame  acceleration  mechanism  is  clearly  "flame  folding"  due  to  the  gradient 
in  the  mean  flow  field  produced  by  the  obstacles.  Iho  tmbulont  burning  velocity 
Sy  is  enhanced  by  an  increase  in  relative  flame  surface  area  (i.o.,  relative 
to  the  projected  flame  area  2»rD)  as  the  flame  is  stretihed  in  the  mean  flow 
gradient  field.  As  the  burning  velocity  increases  the  flow  velocity  ahead  also 
increases  leading  to  more  intense  "flame  folding",  thus  to  a  larger  burning 
velocity  and  so  on.  Superimposed  on  the  large  scale  "flame  fold"  is  the  smaller 
scale  turbulent  flow  field  which  loads  to  flame  wrinkling  and  increased  rate 
of  turbulent,  transport  of  beat  and  mass  across  the  "flume  fold"  surface.  The 
resulting  increase  in  burning  velocity  at  the  "flame  fold"  surface  or  increase 
in  "flame  fold"  burning  rate  leads  to  a  further  increase  in  turbulent  burning 
velocity.  During  these  initial  stages  t ho  flame  acceleration  will  lie  determined 
mainly  by  the  rate  of  "flame  folding".  However,  as  the  flame  accelerates,  the 
smaller  scale  turbulent  flow  field  in  the  wakes  and  shear  layer  regions  will 
begin  to  play  a  more  important  role  in  the  acceleration  process.  As  observed 
in  Figure  b,  the  burning  in  Lite  wake  regions  will  become  less  well  defined  due 
to  smaller  scale  and  larger  intensity  turbulent  entrainment  and  mixing,  and 
the  main  flame  front  will  become  wrinkled  as  the  turbulent p  in  the  wake  regions 
spreads  into  the  main  flow  region. 

Transition  to  turbulence  in  t ho  wake  of  a  cylinder  occurs  at  Reynolds 

numbers  above  200  and  the  location  of  the  transition  point  moves  towards  the 

8  ') 

cylinder  as  the  Reynolds  number  increases.  ’  1  ven  (or  very  low  flow  veloci¬ 

ties  0  m/sec  with  1.2b  cm  obstacles  the  Reynolds  number  is  more  than  HOog 
so  that  the  wake  of  the  obstacles  will  always  be  turbulent.  The  acceleration 
of  the  flame  thus  depends  on  both  the  rate  of  growth  of  the  "flame  fold"  iela- 
tlve  t.o  the  projected  area  2>ir()  and  on  the  increase  in  the  "flame  fold"  burning 
rate.  In  the  initial  stages  of  the  flame  acceleration,  v/here  t  ho  me.u,  flow 
velocity  is  low, the  dominant  mechanism  is  the  rate  of  growth  of  the  "fla-e 
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fold".  Bui  as  the  flame  speed  increases,  the  folded  area  lor  a  cj  i  von  obstacle 
confi  (jura  l  ion  will  approach  some  asymptotic  maxi  mum  since  the  flame  cannot,  he 
stretched  out  indefinitely.  Beyond  this, the  I  lame  acceleration  will  depend 
more  on  the  small  scale  turbulence  intensity  which  qoverns  the  burn i mi  rate 
of  the  "flame  folds".  The  small  scale  turbulence  derives  its  enorqy  from  the  shear 
flow  gradient  of  the  mean  flow  field,  thus  the  flame  speed  will  he  governed 
by  the  rate  in  which  the  fine  scale  turbulence  can  derive  its  energy  from  the 
mean  shear  flow.  For  the  initial  stages  of  the  flame  acceleration  where  the 
development  of  the  "flame  fold"  dominates,  a  computer  simulation  of  the  tran¬ 
sient  inviscid  flow  structure  should  he  adequate  tor  determining  the  flame  ac¬ 
celeration  due  to  "flame  fold"  area  increase.  For  the  later  stages  of  the  flame 
acceleration,  when  the  folded  structure  approaches  some  steady  value  and  small 
scale  turbulence  dominates  subsequent  flame  acceleration,  it  should  he  possible 
to  model  the  acceleration  analytically  by  poslul.ii inq  a  dependence  of  the  burn¬ 
ing  velocity  on  small  scale  turbulence.  In  this  paper  we  shall  propose  a 
simple  feedback  mechanism  for  the  latter  stages  ot  the  flame  acceleration 
based  on  the  assumption  that  the  increase  in  the  rate  ni  burning  of  the  "flame 
fold"  controls  the  acceleration  process.  The  details  of  the  initial  flare 
acceleration  leading  to  i.ho  establ  ishment  of  the  " I  lame  fold"  will  not  he 
considered . 


/l.  R_e_su Its  and  Discussion 

According  to  Townsend'1*  the  turbulent  e  intensify  u  ' /u  in  a  plain' 
wake  is  of  the  order  of  O.d,  where  u1  is  the  root  mean  square  turbulent  velo¬ 
city  and  u()  is  the  mean  velocity  variation  across  the  wale.  if  we  lake  u  as 
one  ha  1 f  the  maximum  flow  velocity  ahead  of  the  flame,  we  obtain  a  turbulent 
velocity  u’  -  0.?  1^.  in  the  wake  of  the  obstacle  ahead  of  the  flame.  Townsend's 
conclusions  are  based  on  fully  developed  turbulent  flow  in  the  far  wake  region 


’  r  'mt*  - * '  ‘vv  a 


r  i  mi  r 


■  of  an  obstacle  in  an  infinite!  stream,  but  the  os  t  imat.e  ol  u '  ()./  II 

ina  x 

(■'•0.2  Rf)  in  the  near  wake  rot) ions  of  Muff  olist.u  It'-,  in  a  boundary  layer  is 

II  )  "> 

supported  by  measurements  by  Chant),  Couni  ban.  Hunt  and  Jackson,  ‘‘  and  by 
recent  measuremonts  in  our  laboratory.  Townsend  also  quote's  a  value  of 
u  ~  O.?l/uo  (i.e.,  0.1  lij  )  for  the  entrainment  velocity,  i.e.,  rate  at 

wliicb  non-turbul out  fluid  flows  through  the  bound intj  surface  of  the  wake  and 
becomes  turbulent.  Thus  as  the  speed  of  the  flame'  increases,  both  the  r.m.s, 
turbulent  velocity  and  the  rate  at  which  the  unhurried  qns  ahead  of  the  flame 
is  made  turbulent  increases. 

To  model  tno  feedback  mechanism  based  on  increase  in  "flame  fold" 
burning  rate  due  to  this  increase  in  turbulence  with  flame  speed,  let  us  focus 
on  the  region  around  the  n^-h  obstacle  and  denote*  the  average  speeds  of  the 
flame  front  through  the  regions  just  before  and  just  after  this  obstacle  by 
R^(n-l)  and  ft,-(n),  respectively.  These  flame  speeds  will  he  determined  by 
the  respective  turbulent  burning  velocities,  which  according  to  our  previous 
discussion  depend  on  both  t.be  relative  si/e  ol  the  "t’.ace  told"  and  on  the 
"flame  fold"  burning  rate.  If  the  relative'  si/e  ol  the  ’  l  lame  fold"  (i.e,, 
relative  to  the  projected  area  ibii'D)  remains  the  same,  the  increase'  in  ilame 
s[)eed  is  due  solely  to  an  increase  in  "flame1  told"  hurtling  rate,  i.e,,  increase 
in  effective  burning  velocity  at  the  "(lame  told'  surtne.e.  Making  l bo  assump¬ 
tion  that  this  burning  rate  is  proportional  to  (o')1'  (tor  the  correlation 
proposed  by  Abdel -Gayed  and  Bradley^  for  high  intensity  isotropic,  turbulence, 
for  example?,  p  =  l).?d!l)  we  find  that 

Rf(n)  =  R f  ( n - 1 )  ( 1  t  r.  ^  '  1  (4) 

where  <s u  ' / u  ‘  is  the  relative  increase  in  turbulence  clue*  to  the  presence  oi  the 
obstacle,  and  we  have  assumed  tliot  u  1  /u  '  1. 


I  be  increase  in  turbulence  in  the  developing  v.aio  region  will  not 


'be  uniform  throughout  the?  region  hetw  on  the  obst  .u:  1  es ,  but.  will  develop  from 
(in  initial  transition  point  wbicb  depends  on  the  Reynolds  numbt'r  of  (In'  flow 
and  tho  initial  turbulence*'’*3  and  spread  as  more  of  Un*  outer  stream  fluid  is 
engulfed  inlo  the  wake,  lor  the  flame  speeds  we  are  concerned  with,  the  ini¬ 
tial  turbulence  and  Reynolds  numbers  are  sul  fie  iently  large  so  tb.it  transition 
can  be  considered  to  occur  right  at  the  obstacle  itself.  Hence,  .su’/u’  which 
represents  the  relative  increase  in  turbulence  averaged  over  the  entire  region 
between  tho  obstacles  will  depend  on  tho  separation  of  the  obstacles  and  on 
the  blockage  ratio  H/0.  If  we  assume  a  linear  growth  of  the  region  with  in¬ 
creased  turbulence  in  tho  near  wake  and  a  power  dependence  on  1  I/P  we  obtain  a 
relation  of  tho  form 


Rf(n)  =  Rf(n-l)!l  <  {',‘f  |  ,  (raJ 

between  tho  average  flame  speeds  before  and  after  the  n *-h  obstacle,  where  u 
is  a  constant  and  L  is  a  length  which  character ,  zes  the  growth  of  the  turbulent 
region.  The  assumption  ol  1  inear  growth  of  the  turbulent,  region  in  tho  near 
wake  region  is  rattier  arbitrary,  but  as  we  shall  see  this  assumption  does  give 
the  correct  scaling  (rom  obstacle  configuration',  with  p  -  3.7i>  cm  to  obstacle 
conf  igurations  wi  th  p  *  10  cm.  A1  tliough  Lquatu.n  ('>)  has  been  obtained  by 
focussing  on  the  regions  just  before  and  just  alter  t  he  i  j  <- » .  obstacle,  the  re¬ 
sult  does  not  depend  on  the  relative  sue  of  tlu>  "tin:-. ■■  told".  As  long  as  the 
relative  size  of  the  "flame  fold"  remains  the  same  and  t ho  relative  increase 
in  turbulence  depends  on  the  obstacle  configuration  only,  exactly  tho  same 
arguments  apply  for  any  size  "flame  fold". 

bincp  the  above  met. ban  i  sin  does  not.  in'  lode  in,.,  (.out  r  i  but  ion  to  flame 
acceleration  due  to  increase  in  the  relative  size  of  the  "flame  fold"  it  can¬ 
not  describe  the  initial  acceleration.  However,  il  we  assume  that  the  int cease 
"flame  fold"  burning  rate  dominates  the  flame  ,v  .  e  1  n  a  l.  i  on  process  alter  (lie 
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(‘lame;  has  propagated  over  n()  obstacle'  ,  then  tor  n  ■  n  w>>  (.an  use  t  •  (ti.-t  1. 1  <  ■ :  >  A 
to  obtain  the  following  relation 


ftf(n)  -  K(  (n())  I  I  i 


,IL  ‘  ("-»„) 
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for  the  flame  velocity  after  the  n^1’  obstacle,  where  ve  have  defined  p  =  i./p. 
For  small  p  the  above  relation  becomes 


(r-r  ) 

II,  (r)  =  k  ( r  )exp[  ■  ° 
t '  I  o  1  l> 


(’Vl 
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(/) 


In  other  words*  it  is  predicted  that  the  t  l.u  velin  it.v  increases  exponentially 
with  r  for  r  >  r  .  As  can  be  seen  in  I  igure  /,  this  prrd  ic  t  i  on  is  supported 

by  the  experimental  results.  If  r  is  chosen  at  GO  <m  with  a  -  0.31  and 

P0  =  37.8  cm.  Equation  7  provides  a  good  correlation  of  the  experimental  results 
from  all  the  experiments  with  obstacle  coni  i  qua  t  ions  which  include  the  !.?!>  cm 
copper  tube  spiral.  In  fact,  for  those  experiments  with  copper  tube  spiral 
obstacles  only,  the  correlation  is  qood  down  to  about,  30  cm  of  flame  travel. 

The  prediction  for  the  larqer  <1  cm  obstacles  with  the  same  parameters  <(  and  p  , 
usituj  l.<|uat.ion  6  rather  than  Iquation  7  since  tin*  exponent  ia  1  approximation 
is  not  valid  in  tliis  cast',  is  shown  in  figure  8.  Aqain,  agreement  with  the 

experimental  results  is  quite  qood  for  r  (it)  tn  (i.e.,  n  '  A),  but  signi.i- 

cant  deviations  are  clearly  observed  lor  r  •  i>)  ii«. 


The  transition  from  flame  acceleration  domino ti>d  by  the  rate  of 
flame  fold  qrowtli  to  acceleration  dominated  by  the  m<  tease  in  "ll.c-e  lolil" 
burn  i  n  <j  rate  is  clearly  not  a  sharp  one.  Their  will  hr  a  t  tans  it.  n  n  reqien 
where  both  mechanisms  are  equally  important,  and  in  tms  reqion  the  t  lame 
acceleration  is  expected  to  be  greatest.  That  the  imti.it  acceleration  of 
the  flame  is  greater  than  that  predicted  by  [.qua  t  ’.mo.  f.  and  1  is  clearly 
evident  in  both  Figures  7  and  H.  Note  that  lor  r  o  all  the  curves  should 
con  ve  i  ge  to  an  initial  flame  v  e  1  oc  i  ty  o  I  3.3  m/  s  .  (in «  e  t.  he  increase  in  "flame 
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.•fold"  burning  rate  l*orj  i m*>  lo  dominate  t.ho  a c*: •* I era t  ion  imin's1,  .1  *.  i tip  1  n  mhh'o  1 
of  the  feedback  mechanism  based  on  a  positive  leedbnk  coupling  between  the 
turbul  once  produced  by  the  obstacles  and  tlie  flame  speed  predicts  an  exponen¬ 
tial  increase  in  flame  speed  which  agrees  well  with  the  rupee bienta 1  resiilts. 

In  fact,  f.quation  '1  which  relates  the  chump;  in  I  I  erne  speed  to  the  champ'  in 
turbulence  can  be  applied  more  generally.  It  would  predict.,  for  example, 
that  the  flame  will  slow  down  if  there  were  no  obstacles  beyond  a  certain  dist¬ 
ance;  the  champ;  in  turbulence  intensity  •Su’  would  then  eventually  become  ne¬ 
gative  and  the  flame  would  decelerate,  furthermore ,  1  hi'  parameters  P(  and 
in  particular  p  ,  would  be  expected  to  depend  on  geometry  (i.o.,  planar,  cy¬ 
lindrical  or  spherical)  and  to  be  constant  only  in  a  certain  range  of  p/ll. 

The  predicted  exponential  increase  in  (lame  speed  cannot  Continue 
indefinitely.  As  the  flame  speed  increases  the  strength  of  the  shock  waves 
produced  will  increase  ami  shock  wave  heating  combined  with  turbulent  mixing 
may  lead  to  auto- ign i  tion  near  the  leading  shod  v.iv,  in  which  case  a  coupled 
shock  wave  reaction  /one  complex  may  continue  t;  pi.mugito  at.  a  <  cud  ant  velo¬ 
city  characteristic  of  the  obstacle  configuration,  'nub  roupled  complexes 
with  velocities  depending  on  the  obstacle  (.out  is  n  a  t  ion  hut.  nurn  valler  than 
the  C-v)  detonation  velocity  have  been  observed  !  .  .■vbuirs.'  '  /•■Hi¬ 

ther  possibility  is  that,  t.ho  rate  0!  chemical  nvi  t  mi  bmomes  the  limitini 
factor.  As  the  flame  speed  increases  the  intensity  ot  the  turbulence  in  tno 
regions  tietween  obstacles  increases,  leading  to  the  observed  exponential  flame 
acccl oration.  However,  there  is  a  maximum  mass  rate  at  which  a  given  voice 
of  a  combustible  gas  at  a  given  pressure  can  bun,  ibis  rorresponds  to  the 
well  stirred  reactor,  ihus  at  a  certain  turbulent  e  intensity  chemistry  lectins 
to  ploy  a  role.  If  tin;  rate  of  turbulent  mixing  exceeds  tin;  rate  at  which 
the?  mixture  can  burn  the  flame  may  be  temporarily  blown  nut  or  locally  guem  bed. 
Such  quenching  has  been  observed  by  Wagner  arid  rol  labor  i  tors  ,  ^ for  exnmole. 


Our  ox  per  i  men  to  1  results  winch  incliil**  I  lame  larger  IJi.m  mOM  m/s.-t  for 

stoichiometric  methane-air  mixtures  show  nn  1  ltd  ica  1. 1  on  t  hut  Lhr  I  I  .imp  speed 
is  approach  i  rig  a  maximum  value.  Thus  wo  would  expect  I  ho  flame  to  rontinue 
to  .accelerate  to  ovpii  larger  flame  speeds  if  t.hr  obst.u  li'  com  i  put  a  Li  on  was 
continued. 


The  maximum  flame  speed  observed  with  no  obstacles  in  the  chamber 
was  rJ  m/sec,  and  no  acceleration  was  observed  in  the  last  i!0  cm  of  propaga¬ 
tion.  This  is  almost  twice  the  value  observed  in  the  to.S  cm  radius  chamber 
in  t ho  laboratory."1^  1  ho  difference  can  tie  attributed  to  additional  surface 
roughness  of  the  wooden  plates  and  slight,  variation  in  plate  separation  due 
to  sagging  of  the  ton  plate.  The  plate  separation  varied  by  as  much  .is  1  cm 
over  the  area  of  the  chamber.  This  variation  of  plate  separation  also  account 
for  the  fact  that  the  flame  speeds  at  30  cm  were  t.\ p ion  1 1  v  larger  than 
those  obseived  in  the  smaller  laboratory  chami-ir  fur  tin*  same  obstacle  i  on- 
figuration.  Tor  11/ 1)  =  tl.3'1  with  1.2b  cm  copper  tubing  obstacles,  for  example, 
the  speed  at.  30  cm  i  •;  -lb  !'  S  nr,  to  he  compare*!  vith  «i l  •  ■  1 1  3‘)  *.i/s  observe.!  in 
the  laboratory  experiments  with  the  same  obstai  les  at  this  li/ll.  The  variation 
in  i>  due  to  sagging  leads  to  an  uncertainty  in  il/l)  of  fl.dn,  and  if  this  un¬ 
certainty  is  taken  into  account  the  two  results  are  in  agreement. 

The  strength  of  the  shock  waves  prodm  ed  by  these  fast,  tuibulent 
flames  were  not  measured  in  a  systematic  manner,  but  t..->  pressure  transducers 
at  the  periphery  of  the  chamber  were  used  to  mm  mil  t  hi*  sheet  wave  overpres¬ 
sures  in  some  of  the  experiments,  i  or  tint  fastest  ilar>>  with  a  maximum  l  lame 
speed  of  'lib  ni/sec  a  shock  wave  overpressure  of  ').(•■!  bar  was  obtained.  Accord¬ 
ing  to  the  steady-stali'  analysis  of  Guirao  el  -s  1  . f>  t  h  i  •.  would  i  on  expend  to 
a  constant  burning  velocity  of  about  3b  m/s  to  he  compared  with  a  mivim.im 
turbulent  burning  velocity  of  bb  m/s  obtained  using  ignation  \  with  a  density 
ratio  of  7.S.  Similarly,  for  the  flame  with  n.ionp  w  lm  ity  .b'-c  -r/s  the 


i'\»  ,  winch  ac<  oi  <!  i  nq  to  human  ot  .1 1  . 


measured  shock  overpressure  was  O.P, 

corresponds  to  ri  constant  hurn  i  n  <  1  velocity  of  ahoul  ,’H  m/s  to  |.o  1  nn.pared 

with  a  ni.i  x  i  mum  turbulent  hurninq  velocity  ol  '11  m/s  oht  .1  i  tied  Icon  iquotion  3. 

In  view  of  t.  ho  dramatic  accel  era  t  ion  ol  tlu'  flames  observed  in  both  of  t  hose 

COSOS  the  cl  i  saqrooment  with  steady  state  calculations  based  on  >  1 ».  1  .im.  1  corn- 

i, 

in<|  velocity  is  not.  sorpi'i  s  i  ruj .  Strehlow  et  n  1  .  hove  recently  shown  that, 
calculations  usinq  constant  velocity  llanos  hosed  rn  tin1  maximum  effective 
hurninq  velocity  of  accelern  tin>i  flames  can  lead  to  an  overestimate  of  the 
maximum  overpressure .  Pur thormore,  Lquation  il  whii.li  we  have  used  to  calculate 
the  turbulent  hurninq  velocity  may  have  to  he  modified  lor  violently  acce 1  era t- 
i  n«l  flames.  1  he  aqreement  with  steady  state  analysis  is  nidi  bet  ter  (or  the 
slower  flames,  i  or  the  flame  with  maximum  speed  01  IVl  m/s,  for  example,, 
the  measured  shock  overpressure  is  O.Uoil  bar  correspond  inq  to  a  steady  state 
burn  ini)  velocity  of  about.  1/  m/s  to  be  compart'd  with  ?t)  m/s  obtained  from 
[‘qua t ion  .1. 

ii.  Conclusion 

1  hi'  influence  of  repeated  ot>st.a<.le\  on  irmly  impandir.q  stou.h 'ewetr i 
nicthane-a  i  r  I  lames  in  cylindrical  qeonetry  ha.  wen  invest,  iqa  ted  in  a  1  lu  'her 
where  t.he  (lame  propaqa  t.  ion  could  he  observed  over  a  distance  of  i.'d  cm.  this 
i  lives  t  i  qa  t  i  on  confirms  the  dramatic  int  1  nonce  el  ■»!■•»!  u  les  on  the  speed  ol 
tlie  flame  reported  previously  and  shows  that  with  rep.-aie.l  obstacles  the  flame 
speed  continues  to  increase  up  t.o  speeds  in  excess  ol  'It id  m/sec.  1  lie  shoe! 
waves  qeneratoif  hy  these  fast  turbulent  t  lames  have  ov>  essures  up  to  d.iel 
bar,  which  could  lead  to  extensive  blast  wave  il.mmqe. 

1  be  mechanism  lesponsihle  for  t.be  rapid  a(  ce  1  .u  a  t  ion  of  the  flame 


in  an  obstacle  environment  can  be  understood  in  terms  ot  the  positive  feed¬ 
back  coup  1  inq  between  the  I  lame  and  the  upstream  t  low  p'ouin.ed  hy  tr.e  specific 


.<  vo  ]  uni'  increase  kmv.s  the  flame.  Hi-  | -hv  s  i  < .»  1  pit  I  u » *  *  nf  (.hi1,  punitive 
feedback  mechanism  proposed  in  a  previous  paper^  is  coni  i  t*n-i*t|  by  spark* 
schlioron  phutopraphs  of  t.  Ik?  I  lame  pn.ipapa  t.  ion  over  obstacles  in  a  im.tan- 
(jul ar  channel. 

Tin;  acu?  1  era  t  ion  of  t.hr:  flan.c?  depends  on  bot.h  the  rate  ot  prov.ih 
of  the  "flame  fold"  relative  to  projected  t  lame  area  and  on  the  increase 
in  the  “flame  fold"  born  i  up  rate.  in  flic'  initial  stapes  nf  the  (lame  ac¬ 
celeration  the  dominant  mechanism  is  the  rate'  ol  pmwth  nf  the  "flame'  told", 
hot  as  the  flame  speed  increases  Lite  rate  of  increase  in  the  "flame  fold' 
burn  i  lip  rate  hepins  to  dominate  the  flame  acceleration  prcu.ess,  A  simple' 
feedback  mechanism  is  proposed  fur  tin?  latter  '.tapi  s  ot  the  flame  acceleration 
based  on  the'  assumption  tli.it  t.la?  increase  in  M rate  nf  t.ari.inp  of  the  "flame 
fold"  due  to  tlie  tin  tail  once  prodoc.  ed  by  tin?  obstacles  controls  the  accelera¬ 
tion  process.  Ibis  model  asstm.es  that  a  "flame  told"  (which  does  not  con¬ 
tinue  to  prow)  has  already  been  established,  ant  that  the  rate  ot  hnrntup  of 
this  "flame  fold"  is  determined  by  the  Uirhim.-i  i  e  field  which  the  I  lame  pio- 
paeja tes  into.  1  ho  details  ot  the  turbulent  re  is  .(.ion  or  Ltn>  nature  of  the 
co.ntiustion  /one:  are  not  considered.  l:.e  mod-1  :  i -’d ;  c.  t  ■.  an  exponential  in¬ 
crease  in  flame  speed,  with  distance  propapa : <> d  .  y.>r  obstacles,  wiik.Ii  is  in 
pood  apreemen  t  with  the  oxper  imentai  results.  Al  t  nnupb  this  exponential  in¬ 
crease  cannot  continue  in  let  i  ni  tolv ,  no  evid.  n,  r  shut  the  flame  speed  is  ap¬ 
proaching  a  limitinp  value  was  tound  in  tin?  pie-.r-nt  experiments. 

fhe  details  ol  t  he  initial  accclcra  t  ...n  ot  the  (lame  which  invol.es 
bo t h  "flame  tolcinuj"  and  increase  in  "flame  fold'  Ournii.p  rate  due  to  turbu¬ 
lence  or  tin'  critical  si/e  of  fhe  "(lame  fold  1  Inr  which  the  rate  oi  >po,-.th 
of  the?  "flame  fold"  can  he  neplected  cannot  be  le  t  .-ih:  iwd  frori  the  present 
analysis.  A  more  detailed  model  ot  the  flow  I  i.'l«l  pre.l-:.  ed  by  1 1»»?  ,\c  ce 1  e?  a  t  - 
inq  flame  propapation  over  obstacles  and  of  the  lie  •'  I  >n  I  u  tone  a  ii.'eraction 


i  is  required  t,o  describe  Llio  initial  flame  accelet  at.  ion  and  t.n  es  1.  im.i  to  LI 
si/e  of  the  flame  fold.  Investigations  with  the  aim  of  (level oping  such  . 
model  are  now  in  progress. 
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Figure  2 


Figure  3 


Figure  4 


Figure  5 


Figure  6 


Schema  Lie  diagram  of  experimental  set-up  For  t  lame 
propaga  t.  i  on  exper  i  ini' n  ts . 

Schema t. i c  cross-sectional  view  of  experimental  chamber 
showing  the  obstacle  conf iguration. 

Schematic  diagram  of  the  flame  propagation  cb.mnol 
with  obstacles. 

Spark-sch 1 i oren  photographs  of  stoichiometric  methane- 
air  flame  propagation  in  a  channel  with  obstacles. 
Obstacle  height  II  =  1.27  cm;  obstacle  separation 
I*  =  3. HI  cm. 

<i )  Channel  width  W  =  0.3  t  (m.  Dry  ice  in  t  tie  second 
and  lourlh  ohtaclos, 

b)  Channel  width  W  -  O.Di’.  cm.  Do  •  1 1  y  ice  used. 

Spark-schl ieron  photographs  oi  stoichiometric  methane- 
air  flame  propagation  in  a  channel  of  width  ?.!>4  cm 
with  obstacles.  Obstacle  height,  li  -  1.2  7  rg  obstacle 
separation  P  -  3.  Hi  cm. 

a)  Dry  ice  in  the  next.  to  last,  nhs'.icle. 

b)  No  dry  ice  used. 

Schematic  diagram  illustrating  t  1  a so  propagation  over 
ohs la c 1 es . 


Figure  7  Name  velocity  Hj  vs.  distance  nl  ,.i  opaga  t  i  on  r  fen' 

stoichiometric  methane-air  Names.  !  lie  solid  curves 
correspond  to  liquation  7  wi  tli  n  ~  0.31  and  p  -  37.3  cm. 
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ci)  Copper  tube  'piral;  II  -  1.2S  cm,  p  =  3./!)  cm  on 

bottom  plate,  plastic,  tube  spiral;  II  =  4  cm  p  =  10  cm 
on  top  plate;  ll/I)  =  ().!>/’. 

b)  Copper  tube  spiral  on  bottom  plate;  ll/l)  =  0.34. 

c)  Copper  tube  spiral  on  bottom  plate;  ll/l)  =  0.2b. 

d)  Copper  tube  spiral  on  bottom  plate;  11/0  =  0.12b. 

The  experimental  point  at  30  cm  is  from  the  laboratory 

results  of  Reference  7. 

Figure  8  Flame  velocity  vs  (iislaiue  of  propagation  r  tor 

stoichiometric  methane-air  (lames.  1  he  solid  curves 
correspond  to  liquation  (>  with  ..  -  0.31  and  p  :  37. b  cm. 
Dashed  curves  are  drawn  tltrouqh  dat.a  points  for  x  -  60  cm, 
a)  Plastic  tube  spirals;  11-4  cm,  p  =  10  cm,  on  both 
top  and  bottom  plates;  11/ b  -  0.(i7. 
h)  Plastic  tube  spiral  on  ho  Mom  plate;  11/0  =  O.b/. 

c)  Plastic  tube  spirals;  ')  ton  plate  only;  t')  top 
and  bottom  plate;  11/0  -  o.4. 

d)  Plastic  tube  spiral  on  not  tom  plate;  l!/0  -  0.33. 
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The  Influence  of  Confinement  on  the  Propagation  of 
Detonations  near  the  Detonahil ity  Limits* 


Paper  submitted  to  the  Liyhteenth  Symposium  (International)  on  Combust 
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1.  Introduction 

One  of  the  fundamental  properties  of  an  explosive  gas  mixture  is 
the  composition  limits  of  detonability..  At  the  present  time  it  is  not  pos¬ 
sible  to  predict  the  detonability  limits  of  a  given  mixture  theoretically 
and  the  limits  must  be  determined  experimentally.  Experiments  to  determine 
the  detonability  limits  are  usually  performed  in  detonation  tubes.  Since  a 
fairly  strong  initiation  source  (eg.,  a  solid  explosive  charge  or  a  volume 
of  readily  detonable  mixture)  is  required  for  near  limit  mixtures,  long  tubes 
are  required  to  ensure  that  the  influence  of  the  initiation  source  on  the 
propagation  of  the  detonation  can  be  neglected  at  the  observation  section 
near  the  end  of  the  tube.  Just  how  long  a  tube  is  required  for  a  given  ini¬ 
tiation  source  is  not  known.  Furthermore,  the  influence  of  the  tube  itself 
(i.e.  ,  the  tube  diameter  and  the  cross-sectional  geometry)  on  the  propagation 
of  the  detonation  is  not  known.  In  fact,  no  generally  accepted  operational 
definition  for  the  detonability  limits  based  on  experiments  in  finite  length 
tubes  has  been  established.  Manson  et  al.^  have  suggested  that  the  fluc¬ 
tuations  of  the  local  detonation  velocity  relative  to  average  velocity 
be  used  as  criteria  for  the  stability  of  the  detonation  in  a  tube.  They  sug¬ 
gest  that  self-sustained  detonations  for  which  U,  and  U  agree  within  about 
i'0.2%  can  be  considered  as  stable.  However,  this  is  an  arbitrary  definition. 

The  appearance  of  near  limit  phenomena  such  as  spin  has  also  been  suggested 

2 

as  a  criteria  for  the  limit  in  a  given  tube,  but  it  has  not  been  established 
that  the  appearance  of  spin  or  other  near  limit  phenomena  in  a  given  tube 
corresponds  to  a  unique  composition. 

Further  studies  are  clearly  required  in  order  to  achieve  a  better 
understanding  of  the  influence  of  the  tube  and  tne  initiation  on  the  near 
limit  propagation  of  detonation  waves  in  tubes  be- fore  a  meaningful  operational 
definition  of  the  limit  can  be  deduced.  This  paper  reports  on  the  first 


results  of  our  investigation  on  near  limit  phenomena  in  tubes  of  different 
diameter. 

2.  General  Considerations 

It  is  now  well  known  that  the  classical  theory  of  a  detonation' 
which  regards  a  detonation  wave  as  a  strictly  one-dimensional  structure  con¬ 
sisting  of  a  shock  wave  followed  by  a  reaction  zone  is  an  inadequate  descrip¬ 
tion  of  the  detonation  phenomena.  Numerous  detailed  investigations  of  the 
structure  of  detonations  over  the  past  50  years  have  shown  that  the  propaga¬ 
tion  of  a  detonation  is  a  complex  three-dimensional  phenomena  involving  the 
interactions  of  finite  amplitude  transverse  waves  with  the  leading  shock  front, 
the  reaction  zone  and  the  boundaries  of  the  system.  The  kinematics  of  these 

interactions  have  been  explored  in  considerable  detail  by  various  investiga- 
3  4  5 

tors.  ’  ’  Although  the  three-dimensional  transverse  wave  structure  of  deto¬ 
nations  is  observed  for  unconfined  detonations,  the  most  detailed  investiga¬ 
tions  of  this  structure  have  been  done  in  confined  rectangular  or  round  deto¬ 
nation  tubes.  In  these  cases,  in  particular  for  conditions  marginal  to  the 
propagation  of  the  detonation  wave  (i.e.,  close  to  the  detonability  limits), 
the  influence  of  the  tube  walls  cannot  be  neglected.  The  tube  walls  have 
two  different  effects;  namely,  an  energy  and  momentum  loss  associated  with 
the  boundary  layers  and  a  stabilizing  effect  on  the  transverse  wave  structure. 
For  small  diameter  tubes  the  observed  decrease  in  velocity  with  decreasing 

tube  diameter*’'  can  be  understood  in  terms  of  the  influence  of  the  boundary 
8  9 

layers.  ’  On  the  other  hand,  it  is  also  observed  that  an  apparently  self- 
sustained  detonation  in  a  confined  tube  fails  once  it  emerges  into  an  area 
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expansion  or  an  unconfined  region, For  a  given  mixture  there  appears 

to  be  a  minimum  critical  tube  diameter  required  in  order  for  the  detonation 

to  continue  to  propagate  under  unconfined  conditions.  It  has  been  suggested 

that  this  critical  tube  diameter  is  related  to  the  characteristic  transverse 

12  13 

wave  spacing  of  the  detonation.  ’  In  other  words,  a  minimum  number  of 
transverse  waves  is  required  for  a  self-sustained  detonation  in  an  unconfined 
situation,  thus  indicating  that  the  pronounced  three-dimensional  structure 
observed  in  tubes  near  the  detonability  limits  is  stabilized  by  the  confine¬ 
ment  provided  by  the  tube  walls. 

A  possible  relation  between  the  spinning  detonations  observed  in 

2 

tubes  and  the  detonability  limits  has  been  discussed  by  Dove  and  Wagner, 
who  suggested  that  the  condition  for  stable  propagation  of  a  detonation  wave 
in  a  tube  is  for  the  reaction  time  to  be  short  enough  to  maintain  the  spin 
mode  of  the  lowest  frequency  in  the  tube.  Thus  if  the  onset  of  the  single¬ 
head  spin  structure  corresponds  to  a  unique  fuel  composition,  the  limit  could 
be  defined  on  this  basis.  Associated  with  this  limit  there  would  then  be 

a  characteristic  chemical  length  scale  which  can  be  related  to  the  tube  dia- 

1 4 

meter  and  geometry  using  the  acoustic  theory  of  spin  detonations  of  Manson 
1 5 

and  Fay.  The  success  of  the  acoustic  theory  in  predicting  the  frequency  or 
pitch  of  the  transverse  or  spinning  vibrations  observed  behind  the  detonation 
front  further  indicates  that  boundary  conditions  do  play  an  important  role 
for  the  propagation  of  detonations  in  confined  tubes.  In  fact,  according  the 
the  acoustic  theory  of  Fay  and  Manson,  the  spin  frequencies  are  entirely  deter 
mined  by  the  boundary  conditions  and  do  not  depend  on  the  details  of  the  coupl 
ing  between  the  gasdynamics  and  the  chemical  kinetics  which  gives  rise  to  the 
transverse  instability  in  the  first  place.  The  only  condition  being  that  the 
reaction  time  or  chemical  time  be  short  enough  to  maintain  the  spinning  mode, 
as  discussed  by  Dove  and  Wagner. 
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The  mechanism  by  which  the  transverse  waves  are  excited  and  main¬ 
tained  is  not  completely  understood.  However,  the  work  by  Barthel  and 
StrehlowJ8  ErpenbeckJ^  Toong^8  and  others  has  clearly  shown  that  acoustic 
and  non-linear  perturbations  can  be  amplified  through  the  coupling  with  che¬ 
mical  energy  release.  The  amplitudes  and  wavelengths  of  the  perturbations 
required  to  trigger  the  various  instablities  are  not  known.  However,  it  ap¬ 
pears  that  detonations  are  unstable  to  perturbations  over  a  fairly  wide  range 

•]  g 

of  wavelengths.  According  to  Toong  the  acoustic  kinetic  interactions  depend 
on  the  order  and  the  enthalpy  of  the  reaction,  the  activation  energy  and  most 
important  of  all  the  ratio  of  the  characteristic  acoustic  time  to  the  chemical 
time.  This  is  further  supported  by  the  work  of  Erpenbeck"' ^  who  established 
the  stability  limits  of  detonations  for  various  degrees  of  overdrive  assuming 
a  first  order  Arrhenius  rate  expression.  Although  the  range  of  wavelengths 
ove*'  which  the  detonation  is  unstable  depends  on  the  activation  energy  and 
the  degree  of  overdrive,  he  finds  that  detonations  are  stable  only  at  short 
wavelengths  (i.e.,  short  compared  to  the  length  of  the  reaction  zone),  and 
also  at  long  wavelengths  for  sufficiently  large  degrees  of  overdrive.  From 
these  investigations  it  appears  that  transverse  waves  with  wavelengths  over 
a  fairly  wide  range  can  be  excited.  Thus  the  transverse  wave  structure  of 
a  detonation  will  depend  on  the  preferred  transverse  mode  (or  modes).  This 
preferred  mode  will  be  determined  nnt  only  by  the  gasdynamic-chemical  kinetic 
coupling,  but  also  by  the  boundary  conditions  (for  example,  the  geometry  and 
diameter  of  the  detonation  tube).  As  long  as  the  characteristic  transverse 
dimensions  associated  with  the  boundary  conditions  are  much  larger  than  the 
characteristic  wavelength  associated  with  the  chemical  kinetics  and  gasdyna- 
mics,  the  boundary  conditions  will  play  a  minor  role  in  determining  the 
transverse  wave  structure.  However,  for  tube  diameters  of  the  order  of  the 
characteristic  transverse  wavelength  or  smaller,  the  boundary  conditions 
will  begin  to  play  a  more  dominant  role,  so  that  for  the  same  mixture  the 
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detonation  phenomena  observed  in  a  small-  diameter  tube  could  be  completely 
different  to  that  which  would  be  observed  in  an  unconfined  situation  or  in 
situations  with  different  boundary  conditions.  Not  only  could  the  structure 
of  the  detonation  wave  be  different,  but  composition  limits  of  detonability 
could  also  vary  with  boundary  conditions.  In  fact,  it  may  be  possible  to 
trigger  "detonation"  phenomena  in  a  tube  outside  of  the  limits  of  detonabi¬ 
lity  for  an  unconfined  situation.  The  phenomena  of  "galloping"  detonations 

19  20  21 

observed  by  Mooridan  and  Gordon,  Manson  et  al,  and  Edwards  and  Morgan 

may  be  an  example  of  such  a  phenomena.  If  this  is  the  case  and  if  the  onset 
of  the  "galloping"  mode  is  sufficiently  precise,  then  the  onset  of  the 
"galloping"  mode  could  also  provide  a  criteria  for  determining  the  detonabi¬ 
lity  limits.  The  "galloping"  mode  is  a  longitudinal  mode  with  periodic  des¬ 
truction  and  reformation  of  the  detonation  and, as  has  been  pointed  out  by 

22 

Urticw  and  Oppenheim,  the  reformation  process  is  identical  to  the  transi¬ 
tion  form  deflagration  to  detonation.  "Galloping"  detonations  can  therefore 
be  considered  to  consist  of  periodic  transitions,  in  which  case  the  tube 
walls  and  confinement  are  known  to  play  an  important  role.  However,  the  role 
of  the  transverse  waves,  which  are  observed  in  near  limit  mixtures  in  main¬ 
taining  the  detonation  wave  is  not  understood. 

We  have  undertaken  a  program  to  investigate  the  propagation  of 
detonations  in  long  tubes  (14  m)  of  different  diameters  (5  cm  up  to  30  cm). 
The  aim  of  the  program  is  to  determine  a  criteria  for  establishing  the  deto¬ 
nability  limit  in  unconfined  situations  based  on  laboratory  experiments  and 
to  investigate  near  limit  phenomena  to  clarify  the  role  of  confinement  on 
the  propagation  of  detonations.  This  paper  reports  on  the  first  part  of 
this  investigation  which  involved  a  detailed  study  of  the  propagation 
of  detonations  in  ethylene-air  mixtures  near  the  lean  limit  in 
two  long  tubes  (14  m)  of  diameter  5  and  13  cm.  The  transverse  wave  structure 
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of  the  detonations  was  observed  by  pressure  transducers  at  various  positions 
along  the  tubes  as  described  in  detail  in  the  next  section.  A  detailed  dis¬ 
cussion  and  interpretation  of  the  results  are  given  in  Sec.  4  and  Sec.. 5  is 
the  conclusion. 

3.  Experimental  Details 

The  experimental  apparatus  consisted  of  two  steel  tubes  of  4.8  cm 
and  14.5  cm  in  diameter  with  respective  lengths  of  14.15  m  and  14.63  m.  The 
ethylene-air  gas  mixture  was  prepared  in  continuous  flow  with  the  required 
flow  rates  of  each  component  gas  monitored  via  the  pressure  drop  across  ca¬ 
pillary  tubes. 

The  accuracy  of  the  gas  mixture  produced  was  within  0.1%  for  the 
flow  system  associated  with  the  4.8  cm  diameter  tube  and  within  0.2%  for 
the  14.5  cm  tube.  To  ensure  proper  ga-s  concentration  in  the  experimental 
apparatus,  the  filling  procedure  consisted  in  first  evacuating  the  tubes, 
then  filling  them  with  one  atmosphere  of  the  required  gas  mixture,  followed 
by  a  further  flowing  of  gas  mixture  through  the  tubes  at  one  atmosphere  for 
an  equivalent  of  five  fill  times.  All  the  experiments  were  done  at  one  at¬ 
mosphere. 


The  velocity  and  pressure  records  were  obtained  from  P.C.B. 
Piezotronics  transducers  located  at  different  positions  along  the  tubes 
(Fig.  1).  In  the  4.8  cm  tube  the  first  two  transducers  were  located  at 
2.3  m  and  3.3  m  from  the  ignition  flange.  Next  a  series  rr  five  transdu¬ 
cers  0.5  m  apart,  with  the  first  being  used  as  a  trigger,  was  located  near 
the  middle  of  the  tube  starting  at  a  distance  of  7.32m  from  the  ignition. 

Six  more  transducers  were  located  near  the  end  of  the  tube.  The  first  trans¬ 
ducer  situated  at  10.82  m  was  u^ed  as  a  trigger.  The  next  transducer  was 
connected  to  a  Biomation  digital  recorder  so  that  clear  pressure  traces  could 
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be  obtained.  All  other  transducers  were  used  in  conjunction  with  standard 
oscilloscopes  to  obtain  their  respective  pressure  traces.  Similarly  for 
the  14.5  cm  tube  a  total  of  nine  transducers  were  used,  with  three  located 
near  the  middle  starting  at  a  distance  of  7.32  m  and  separated  by  a  distance 
of  60.96  cm.  The  final  six  transducers  were  located  near  the  end  of  the' 
tube  as  indicated  in  Figure  lb. 

The  ignition  system  consisted  of  flowing  a  slug  of  oxy-acetylene 
gas  mixture  into  the  tube  and  igniting  it  with  a  high  voltage  capacitor 
spark.  The  gas  mixture  concentration  was  varied  between  equimolar  and  20% 

-  80%  Og,  and  slug  lengths  between  three  meters  and  half  meter.  A  slug 
of  half  a  meter  was  found  to  be  the  minimum  length  in  which  direct  initiation 
of  a  detonation  could  be  reliably  obtained  because  of  diffusion  effects  bet¬ 
ween  the  oxy-acetylene  and  the  test  gas  mixture  already  in  the  tube.  The 
characteristics  of  the  ignitor  were  determined  by  detonating  different  con¬ 
centrations  and  lengths  of  the  ignitor  slug  into  air.  The  shock  velocities 
were  measured  at  the  end  of  the  tube  for  the  different  ignitors  and  in  all 
cases  relatively  weak  shocks  of  similar  strength  (450-550  m/s)  were  observed. 

Finally  two  spirals  to  be  used  in  the  2.85  cm  diameter  tube 
were  made  from  0.32  cm  copper  tubing.  The  length  of  each  was  1.5  m  and  the 
pitches  were  15  cm  and  1.6  cm.  They  could  be  positioned  either  near  the 
beginning  (1.535  m)  or  near  the  middle  (6. 01). of  the  tube. 

4.  Results  and  Discussion 

The  detonation  velocities  for  different  compositions  of  ethylene- 
air  observed  near  the  end  of  the  tubes  are  compared  with  the  theoretical 
C-J  value  in  Fig.  2.  The  experimental  points  correspond  to  the  average 
velocity  measured  from  the  set  of  transducers  near  the  end  of  tubes.  As 
can  be  seen  from  Fig.  2,  the  velocities  measured  in  the  4.8  cm  diameter 
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tube  are  in  good  agreement  with  the  theoretical  values  over  the  whole  range 
of  compositions  from  6%  to  3%  ethylene  in  CgH^-air.  However,  in  the  larger 
14.5  cm  diameter  tube  the  observed  velocities  begin  to  consistently  deviate 
significantly  from  the  theoretical  value  below  5%  in  C2H^-air.  At 
3.5%,  for  example,  the  observed  average  velocity  is  approximately  1 3"i  lower 
than  the  theoretical  value.  As  will  be  discussed  later,  this  region  also 
corresponds  to  the  region  where  the  pressure  structure  of  the  detonations 
in  the  two  tubes  differs,  significantly. 

These  results  are  for  "self-sustained"  detonations  which  are 
directly  initiated  by  the  detonation  of  the  acetylene-oxygen  charge  (super¬ 


critical  initiation)  and  whose  observed  veloo't}  is  constant  over  the  last 


half  of  the  tubes  within  experimental  errors  (t  3%).  If  the  initiation 
charge  is  not  strong  enough  to  produce  a  detonation  initially  (i.e.,  sub- 
critical  initiation)  so  that  a  decoupled  shock  flame  complex  is  observed  at 
the  first  two  transducers  in  the  4.8  cm  diameter  tube,  or  if  obstacles  in 
the  form  of  Shchelkin  spirals  are  placed  in  the  4.8  cm  diameter  tube,  com¬ 
pletely  different  phenomena  are  observed.  In  these  cases  large  fluctuations 
in  velocities  are  observed  for  mixtures  which  exhibit  strong  single  head 
spin  structure  for  super-critical  initiation.  The  phenomena  is  analogous  to 
the  "galloping"  detonations  observed  by  Mooridan  and  Gordon,^  Manson  et  al 
and  Edwards  and  Morgan.  The  wave  appears  to  propagate  in  a  cyclic  manner  i 

exhibiting  large  velocity  fluctuations  with  velocities  ranging  from  2100  m/s 
down  to  about  900  m/s.  This  phenomena  of  cyclic  propagation  is  observed  for 

•I 

subcritical  initiation  for  3.5%  02^  and  below.  Unfortunately,  with  our  pre¬ 
sent  ignitor  system,  subcritical  initiation  above  3.5%  C2H4  could  not  be 
achieved  in  a  reliable  manner,  so  that  it  cannot  be  concluded  that  similar 
propagation  phenomena  are  not  possible  in  richer  mixtures.  When  Shchelkin 
spirals  placed  near  the  beginning  of  the  tube  were  used  to  dissociate  the 
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detonation  wave,  the  observed  transition  near  the  end  of  the  tube  is  the 
same  as  that  of  a  galloping  cycle,  and  for  fuel  compositions  below  about 
4.5%  the  original  spinning  detonation  is  not  re-established  by  the  end 
of  the  tube.  Thus,  aithough  apparently  stable  detonations  are  observed  in 
the  4.8  cm  diameter  tube  all  the  way  down  to  3 %  in  C2H^-air,  the  de¬ 
tonations  which  exhibit  clear  single  head  spin  are  very  sensitive  to  both 
the  initial  conditions  and  to  the  boundary  conditions.  By  changing  the  ini¬ 
tiation  energy  or  by  placing  obstacles  in  the  tube  completely  different  modes 
of  propagation  are  observed.  Also,  detonations  in  the  same  near-limit  fuel 
mixture  but  in  the  two  different  diameter  tubes  have  significantly  different 
structure,  c  ■  the  propagation  velocity  appears  to  be  constant  in 

both  tubes. 


For  mixture  compositions  far  removed  from  the  detonability  limit 
the  pressure  records  show  a  typical  multi -headed  detonation  in  both  the  small 
(4.8  cm  diameter)  and  the  large  (14.5  cm  diameter)  tube.  Examples  of  such 
pressure  records  are  shown  in  Fig.  3.  Even  for  these  compositions,  chore  are 
distinct  pressure  variations  of  the  order  of  5-10  atm  near  the  front  result¬ 
ing  from  the  three  dimensional  structure  of  the  detonations.  However,  the 
frequencies  of  these  variations  (MOO  Kc/s)  are  an  order  of  magnitude  larger 
than  the  characteristic  lowest  mode  spin  frequency  (MO  Kc/s)  and  the  pres¬ 
sure  vibrations  disappear  within  two  or  three  cycles.  Thus  for  these  mix¬ 
tures  the  characters sti c  chemical  gasdynamic  transverse  wavelength  is  an 
order  of  magnitude  smaller  than  the  dimensions  of  the  tube  and  there  is  ap¬ 
parently  no  strong  coupling  of  the  detonation  to  the  boundary  conditions. 

This  should  be  contrasted  with  the  structure  observed  for  near-limit  mixtures. 
Fig.  4  shows  typical  examples  of  pressure  records  observed  for  3.3% 

Notice  that  the  structure  of  the  detonation  in  the  small  tube  (Fig.  4a)  dif¬ 
fers  markedly  frcm  that  observed  in  the  large  tube  (Fig.  4b).  Both  pressure 
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traces  show  the  characteristic  pressure  variations  due  to  a  single-head  spin¬ 
ning  wave  behind  the  detonation  front,  but  since  the  diameter  of  the  tubes 
are  different  the  spin  frequencies  are  different.  In  the  small  tube  the 
characteristic  time  of  the  pressure  oscillations  behind  the  front  ot  =  0.1  ms 
corresponding  to  spin  pitch  P  =  6tU  (where  U  is  the  wave  velocity)  of  about 
15  cm  and  a  spin  pitch  to  tube  diameter  ratio  P/D  =  3.1.  Similarly,  in  the 
large  tube  the  characteristic  time  St  =  0.27  ms  corresponding  to  a  P/D  ratio 
of  about  2.5.  The  pressure  pulses  observed  in  Fig.  4  correspond  to  pressure 
variations  of  the  order  of  10  atm  and  are  clearly  not  acoustic  pulses,  never¬ 
theless  it  is  useful  to  compare  the  observed  ratio  of  P/D  to  that  obtained  by 
considering  the  pulses  to  be  due  to  an  acoustic  spinning  wave.  For  a  purely 

transverse  acoustic  wave  we  have^  P/D  =  ?—  (-),  where  c  is  local  speed  of 

Kn  c 

sound  and  k  is  the  root  of  the  first  derivative  of  the  Bessel  function  of 

— — r 

order  n  (k-[  =  1.841  and  kg  =  3.054).  The  theoretical  value  at  3.3'X  CpH^  is 
P/D  =  2.97  which  is  in  good  agreement  with  the  pitch  to  diameter  ratios  ob¬ 
served,  especially  in  view  of  the  approximations  involved  (i.e.,  acoustic 
waves,  C-0  values  for  U  and  c).  From  these  observations  it  is  clear  that  for 
this  composition  the  spin  pitch  or  spin  frequency  is  determined  by  the  cha¬ 
racteristic  dimensions  of  the  tubes  rather  than  by  some  characteristic  length 
associated  with  the  rate  of  chemical  reaction.  This  conclusion  is  further 
supported  by  the  observations  that  in  the  small  tube  the  same  spin  pitch  to 
diameter  ratio  is  observed  over  the  range  of  compositions  from  4.2"  to  3" 
CgH^,  indicating  that  for  these  compositions  the  confinement  or  boundary  con¬ 
ditions  determine  the  transverse  modes  which  dominate  the  propagation  of  the 
detonation.  Similarly,  in  the  large  tube  transverse  modes  characteristic  of 
the  tube  rather  than  the  mixture  dominate  the  propagation  between  about  3.5 
and  3«  CgH^.  At  3"  C2H4  some  of  the  single  head  spin  detonations  that  ap¬ 
peared  to  be  stable  over  the  first  half  of  the  tubes  (i.e.,  for  about  7  m) 
were  observed  to  decay  completely  by  the  end  of  the  cube.  Thus  as  the 
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characteristic  spinning  mode  associated  with  the  tube  gets  further  and  fur¬ 
ther  out  of  tune  with  the  characteristic  chemical  time,  or  equivalently  as 

2 

stated  by  Dove  and  Wagner,  the  chemical  reaction  time  gets  too  long  to 
maintain  the  spin  mode,  the  detonation  becomes  unstable.  Below  3«  we 
were  unable  to  maintain  any  form  of  detonation  modes  in  either  of  the  tubes. 

As  mentioned  previously  the  single  head  spin  mode  of  propagation 
is  also  very  sensitive  to  the  initial  conditions,  indicating  that  the  coupl¬ 
ing  between  the  gasdynamics,  the  chemical  energy  release  and  the  boundary 
conditions  must  be  properly  established  for  a  steady  state  detonation  to 
propagate  the  length  of  the  tube.  If  this  coupling  is  not  established  from 
the  beginning  a  "galloping"  mode  of  propagation  is  observed.  This  phenomena 
appears  to  be  controlled  by  the  amplification  of  transverse  waves  in  the 
reaction  zone  behind  a  leading  shock  front.  Pressure  records  of  the  dif¬ 
ferent  phases  of  the  propagation  for  3*3%  C2H4  in  the  4.8  cm  diameter  tube 
are  shown  in  Fig.  5.  Fig.  5a  shows  a  fairly  weak  shock  wave  (Ap  =  6  atm) 
followed  by  a  reaction  zone  with  intense  pressure  oscillations  in  excess  of 
30  atm  peak  to  peak  with  a  frequency  of  the  order  of  the  spin  frequency.  The 
second  trace  which  is  taken  50  cm  downstream  of  the  first  trace  shows  the 
amplification  of  the  pressure  oscillations  and  the  catching-up  of  the  reac¬ 
tion  zone  to  the  leading  shock  wave.  The  velocity  of  the  leading  shock  wave 
is  966  m/s  or  about  63%  of  the  theoretical  C-J  velocity.  The  stage  at  which 
the  transverse  oscillations  have  just  caught  up  to  the  leading  front  is  shown 
in  Fig.  5b.  This  leads  to  an  overdriven  detonation  with  a  velocity  of 
2087  m/s  which  Uten  decays  to  1687  m/s  between  the  last  two  pressure  trans¬ 
ducers.  (Both  the  two  top  and  two  bottom  traces  are  taken  50  cm  apart,  and 
the  two  middle  traces  are  taken  1  m  apart).  A  pressure  trace  taken  50  cm 
prior  to  the  first  trace  of  this  figure  shows  a  weak  shock  front  (■=  7  atm) 
followed  by  intense  pressure  variations  at  a  time  interval  of  0.25  ms  behind 


the  leading  front.  Notice  that  the  transverse  wave  is  absent  in  the  over¬ 
driven  state.  It  is  not  clear  why  it  has  dissapeared  so  quickly,  however 
the  effective  reaction  zone  would  now  be  much  smaller  and  the  low  frequency 
mode  would  be  incapable  of  being  further  exited.  This  is  supported  by 
Erpenbeck's^  calculation  showing  that  overdriven  detonations  are  stable 
to  a  wider  range  of  wavelengths  than  C-J  detonations. 

The  overdriven  or  fast  mode  of  propagation  is  further  illustrated 
in  Fig.  6a.  In  this  case  the  velocity  between  the  two  top  traces  (separated 
by  50  cm)  is  2050  m/s,  and  the  wave  decays  to  1650  m/s  between  the  last  two 
traces  (separated  by  1.5  m).  Notice  that  the  characteristic  pressure  oscil¬ 
lations  begin  to  reappear  as  the  velocity  approaches  the  C-J  value  of  about 
10  m/s,  lending  support  to  our  previous  argument  for  the  absence  of  trans¬ 
verse  wave  in  the  overdriven  state  based  on  Erpenbeck ' s^ 7  analysis.  The 
pressure  trace  shown  in  Fig.  6b  is  taken  at  a  position  between  the  bottom 
two  traces  of  Fig.  6a.  This  trace  shows  that  the  wave  actually  consists  of 
a  relatively  strong  shock  front  (&p  =  11.2  atm)  followed  by  a  high  pressure 
region  (ap  -  20  atm) . 

These  results  clearly  show  that  transverse  waves  which  are  coupled 
to  the  tube  are  amplified  in  the  chemical  reaction  zone  behind  the  leading 
shock  wave  and  a  mode  of  unsteady  propagation,  consisting  of  transverse  wave 
amplification  with  subsequent  catch-up  of  the  transverse  wave  region  to  the 
shock,  resulting  in  an  overdriven  detonation  which  then  decays  and  leads  to 
a  separation  of  the  shock  wave  and  the  reaction  zone.  The  process  then  re¬ 
peats  itself.  Transverse  wave  structure  in  the  shock  flame  region  during 

the  dissociated  or  slow  phase  of  galloping  were  also  observed  by  Manson  et 
20 

al ,  but  their  role  in  the  unsteady  propagation  was  not  recognized.  The 
explosion  which  occurs  as  the  transverse  wave  catches  up  to  the  leading 
shock  wave  is  similar  to  that  which  can  also  occur  in  flame  to  detonation 
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transition  as  Urtiew  and  Oppenheim  have  pointed  out. 

To  further  clarify  the  role  of  the  transverse  waves  for  the  propa¬ 
gation  of  detonations  in  tubes  we  inserted  Shchelkin  spirals  at  various  po¬ 
sitions  in  the  small  tube.  Two  spirals  of  different  pitch  were  used,  one 
of  pitch  15  cm  which  corresponds  approximately  to  the  spin  pitch  observed 
in  this  tube.  It  was  found  that  this  spiral  had  almost  no  influence  on  the 
propagation  observed  at  the  end  of  the  tube.  However,  for  the  spiral  of 
pitch  1.6  cm  the  effect  is  quite  dramatic.  As  shown  in  Fig.  7a  this  spiral 
completely  destroys  the  transverse  wave  structure  and  the  spinning  detonation 
is  completely  dissociated  as  it  leaves  the  spiral  region  as  shown  on  the 
bottom  trace.  When  this  spiral  is  placed  near  the  beginning  of  the  tube, 
reformation  due  to  the  build-up  of  transverse  waves  occurs  for  composition 
above  3.5%  C2H4.  A  typical  record  shows  the  transition  about  7  m  down¬ 
stream  of  the  spiral  regicn  for  4%  C2H^  is  shown  in  Fig.  7b.  This  transi¬ 
tion  is  virtually  identical  to  one  phase  of  the  galloping  propagation  ob¬ 
served  with  sub-critical  initiation  in  3.3%  C^ll^  shown  in  Fig.  7c.  Notice 
that  in  both  cases  (i.e..  Fig.  7b  and  7c)  the  transverse  vibrations  in  the 
reaction  zone  amplify  and  finally  merge  with  the  leading  shock  wave.  In 
both  cases,  higher  frequency  modes  are  excited  behind  the  leading  shock 
wave  and  the  velocities  of  the  leading  shock  waves  are  between  900  and  1200  m/s 
(i.e.,  less  than  the  C-J  detonation  velocity)  during  the  merging  process. 

The  near  limit  detonation  phenomena  observed  in  confined  tubes 
are  clearly  strongly  influenced  by  the  spin  vibrations  associated  with  the 
characteristic  modes  of  the  tube.  In  fact,  for  mixtures  whose  characteris¬ 
tic  gasdynamic-chemical  transverse  wavelength  is  less  than  or  of  the  order 
of  the  transverse  dimensions  of  the  tube, completely  different  modes  of  pro¬ 
pagation  are  observed  depending  on  the  coupling  between  the  gasdynamics, 
chemical  energy  release  and  the  walls  of  the  tube.  Due  to  the  strong 
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influence  of  the  confinement,  both  single  head  spin  and  galloping  detona¬ 
tions  are  observed  over  a  fairly  wide  range  of  mixture  compositions.  These 
near  limit  phenomena  depend  more  on  the  tube  than  on  the  chemical  kinetics 
of  the  mixture.  For  different  size  or  geometry  tubes  these  near  limit 
phenomena  would  therefore  be  expected  to  occur  at  different  mixture  compo¬ 
sitions.  Thus  detonability  limits  which  are  to  be  extrapolated  to  infinite 
tube  diameter  or  to  the  unconfined  situation  must  exclude  these  phenomena. 

We  therefore  propose  that  these  detonability  limits  be  characterized  by  the 
composition  at  which  clear  single-head  spin  first  appears  in  the  tube. 

At  this  composition  the  characteristic  gasdynamic-chemical  transverse  wave¬ 
length  is  approximately  tuned  to  the  size  of  the  tube.  This  does  not  answer 
the  question  of  the  number  of  transverse  waves  required  for  an  unconfined 
detonation,  but  it  does  provide  a  length  scale  for  a  given  mixture  composi¬ 
tion,  i.e.,  the  tube  diameter  or  the  pitch  of  the  spin,  which  can  be  used 
to  estimate  the  characteristic  wavelength  in  an  unconfined  situation  for 
this  mixture  composition. 

If  we  use  the  above  criterion,  we  find  that  the  lean  detonability 
limit  for  ethylene-air  mixtures  lies  between  4.5  and  4.2%  in  the  4.8  cm 
diameter  tube  and  between  4  and  3.5%  C2H4  in  the  14.5  cm  diameter  tube.  The  cha¬ 
racteristic  transverse  wavelength  of  the  mixture  has  therefore  increased  by  a  factor  of 
about  3  for  a  composi tional  change  of  about  0.7%  C^H^.  Typical  pressure  re¬ 
cords  showing  the  appearance  of  single  head  spin  between  these  compositions 
are  shown  in  Fig.  8.  Fig.  8a  and  b  are  pressure  records  of  detonations  in 
the  4.8  cm  diameter  tube  at  4. 5'%  and  4.2%  CgH^,  respectively  (note  the  dif¬ 
ference  in  vertical  scale).  At  4.5%  the  characteristic  time  of  the  pressure 
oscillations  is  of  the  order  of  50  ysec  (with  still  higher  frequencies  ob¬ 
served  near  the  front).  A  time  of  50  nsec  gives  a  spin  pitch  of  about  8.35 
cm  and  a  pitch  to  diameter  ratio  of  about  1.72.  This  therefore  corresponds 
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to  the  second  lowest  spin  mode  whose  pitch  to  diameter  ratio  calculated  from 
acoustic  theory  is  1.81.  At  4.2 «,  the  character!- Stic  time  is  about  90  ns 
(although  there  is  some  slight  variation  in  frequency  near  the  front)  giving 
a  pitch  to  diameter  ratio  of  about  3,  which  corresponds  to  the  lowest  spin 
mode  in  this  tube.  Fig.  8c  and  d  are  the  pressure  records  at  4%  and  3.5ii 

in  the  14.5  cm  tube,  respectively  (note  the  different  in  vertical  scale). 
The  pressure  oscillations  near  the  front  at  4»  are  separated  by  a  time  of 
about  0.1  ms  to  0.15  ms,  which  is  about  half  of  the  characteristic  time  as¬ 
sociated  with  the  lowest  spin  mode  in  this  tube,  indicating  that  a  higher 
mode  has  been  excited,  at  least  near  the  front.  This  can  be  compared  with 
the  pressure  oscillations  at  2.5%  which  show  no  indication  of  any  mode  other 
than  the  lowest  spin  mode  whose  characteristic  time  is  about  0.28  ms,  cor¬ 
responding  to  a  pitch  to  diameter  ratio  of  about  2.8. 

The  variation  in  characteristic  wavelength  with  composition  is 
expected  to  be  very  dramatic  near  the  detonability  limits  and  it  is  expected 
that  for  sufficiently  large  diameter  tubes,  single  head  spin  would  not  be 
possible.  The  critical  tube  diameter  above  which  single  head  spin  is  not 
observed  would  then  correspond  to  maximum  transverse  wavelength  possible 
for  the  mixture,  and  the  composition  corresponding  to  this  characteristic 
wavelength  would  be  the  composition  limit  beyond  which  a  detonation  cannot 
be  initiated  no  matter  how  large  a  volume. of  mixture  is  available.  Clearly 
experimental  investigations  in  different  diameter  tubes  are  required  before 
any  extrapolation  can  be  attempted.  These  investigations  are  now  in  progress. 

5.  Concl usion 

The  propagation  of  detonations  in  confined  tubes  is  clearly  in¬ 
fluenced  by  the  confinement  provided  by  the  tube  walls.  This  influence 
becomes  particularly  dramatic  for  marginal  detonation  waves  with  characteristic 
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transverse  wavelengths  of. the  order  of  or  less  than  the  tube  diameter.  In 
these  cases  single  head  spinning  detonations  whose  pitch  is  determined  by 
the  diameter  of  the  tube  are  observed  over  a  fairly  wide  range  of  composi¬ 
tions.  The  structure  of  detonation  therefore  depends  critically  on  the  tube 
diameter  and  is  maintained  by  the  interaction  of  a  combustion  zone  instabi- 

♦ 

lity  which  couples  to  the  characteristic  modes  of  the  tube.  Thus  for  dif¬ 
ferent  diameter  tubes,  different  characteristic  transverse  wavelengths  are 
observed  for  the  same  mixture  composition.  Since  these  spinning  modes, 
which  are  strongly  coupled  to  the  tube  dimensions,  would  not  occur  at  the 
same  composition  under  unconfined  conditions  or  in  tubes  of  larger  diameter 
we  have  proposed  that  the  equivalent  unconfined  detonability  limits  be  cha¬ 
racterized  by  the  composition  at  which  single  head  spin  is  first  observed  in  the  tube. 
This  then  provides  a  length  scale  associated  with  the  detonability  limit 
which  can  be  extrapolated  to  different  situations.  Based  on  this  criteria 
we  obtain  a  detonability  limit  between  4. 5-4. 2%  ethylene  in  CgH^-air  in  a 
4.8  cm  tube  and  a  limit  between  4-3.5%  in  a  14.8  cm  tube.  Although  this  is 
clearly  insufficient  for  any  extrapolation  it  does  provide  a  basis  for  fur¬ 
ther  work  in  different  size  tubes  which  is  now  in  progress. 

The  role  of  transverse  spinning  waves  in  maintaining  a  detonation¬ 
like  phenomena  in  a  round  tube  have  been  further  clarified  by  varying  the 
initiation  energy  and  by  inserting  Shchelkin  spirals  at  various  positions 
along  the  4.8  cm  diameter  tube.  It  is  found  that  waves  which  exhibit  appa¬ 
rently  stable  single  head  for  the  length  of  the  tube  for  super-critical  ini¬ 
tiation  propagate  in  an  unsteady  cyclic  fashion  analogous  to  galloping  de¬ 
tonations  for  sub-critical  initiation.  The  propagation  of  these  unsteady 
waves  is  controlled  by  the  amplification  of  transverse  vibrations  in  the 
reacting  zone  behind  the  leading  shock  wave.  As  the  transverse  vibrations 
amplify  they  catch  up  to  the  shock  wave  leading  to  an  overdriven  detonation. 
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In  this  overdrive  stage  the  transverse  vibrations  disappear  and  only  begin 
to  appear  again  as  the  velocity  approaches  the  C-J  value.  The 
process  then  repeats  itself.  A  similar  phenomena  is  observed  with  the 
Shchelkin  spiral  in  the  tube.  In  this  case  the  spiral  destroys  the  trans¬ 
verse  wave  structure  and  the  wave  dissociates,  transverse  waves  then  amplify 
in  the  reaction  zone  region  and  transition  occurs  when  these  waves  catch 
up  to  the  leading  flame  front. 

At  this  stage  of  our  investigation  the  results  are  qualitative  in 
nature  but  our  detailed  observations  on  the  role  of  transverse  vibrations  in 
maintaining  detonation-like  phenomena  in  confined  tubes  suggests  a  new  ap¬ 
proach  for  further  experimental  and  theoretical  investigations.  These  inves¬ 
tigations  are  now  in  progress. 
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Fig.  1  Schematic  diagram  of  the  experimental  apparatus. 

a)  4.8  cm  diameter  tube, 

b)  14.5  cm  diameter  tube. 


Fig.  2  Comparison  of  average  detonation  velocities  measured  at 

the  end  transducers  in  the  two  detonation  tubes  with 
theoretical  Chapman-Jouguet  detonation  velocities  for 
different  fuel  compositions. 

Fig.  3  Typical  pressure  records  of  multi-head  detonations  in 

ethylene-air  mixtures. 

a)  5%  in  4.8  cm  diameter  tube.  Biomation  record: 
Time  base:  20  usec/div. 

Vertical  scale:  6.8  atm/div. 

b)  7.8%  in  14.5  cm  diameter  tube.  Biomation  record: 
Time  base:  20  usec/div. 

Vertical  scale:  6.8  atm/div. 


Fig.  4  Pressure  records  showing  single-head  spin  detonations  in 

ethylene-air  mixtures. 

a)  3.3%  02^^  in  4.8  cm  diameter  tube.  Biomation  record: 
Time  base:  100  us/div. 

Vertical  scale:  6.8  atm/div. 

b)  3.3%  in  4.8  cm  diameter  tube.  Biomation  record: 
Time  base:  100  us/div. 

Vertical  scale:  6.8  atm/div. 

Fig.  5  Pressure  records  of  near  limit  phenomena  observed  in  the 

4.8  cm  diameter  tube  for  sub-critical  initiation  of  3.3% 
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a)  Amplification  of  the  spinning  wave  behind  the  shock 
front  at  the  end  transducers. 

Time  base:  0.2  ms/div. 

Vertical  scale:  13.6  atm/div. 

b)  The  transition  process  observed  at  the  end  transducers. 

Time  base:  0.2  ms/div. 

Vertical  scale:  13.6  atm/div. 

Pressure  records  showing  the  fast  phase  of  the  propagation 
of  galloping  detonation  for  sub-critical  initiation  at 
3.3  C2H4  in  the  4.8  cm  diameter  tube. 

a)  Top  three  traces  are  taken  0.5  m  apart  and  the  bottom 
two  traces  are  taken  1.5  m  apart  (at  the  end  transducers) 
Time  base:  0.5  ms/div. 

Vertical  scale:  13.6  atm/div. 

b)  Biomation  record  taken  between  two  bottom  traces  in  Fig.  6a. 
Time  base:  20  ysec/div. 

Vertical  scale:  3.4  atm/div. 

a)  Pressure  records  taken  in  the  region  of  spiral  obstacles 
of  pitch  1.6  cm  in  the  middle  of  the  4.8  cm  diameter  tube 
showing  the  dissociation  of  the  spinning  wave  by  the 
obstacles  at  3.5%  C^. 

Time  base:  0.5  ms/div. 

Vertical  scale:  11.3  atm/div. 

'  b)  Pressure  records  showing  the  transition  process  7m 
downstream  of  spiral  obstacles  of  pitch  1.6  cm  for  4% 

C2H4  in  4.8  cm  diameter  tube. 

Time  base:  0.2  ms/div. 

Vertical  scale:  11.3  atm/div. 
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c)  Pressure  records  taken  at  the  middle  transducers  for 
sub-critical  initiation  of  3.3%  in  the  4.8  cm  tube. 
Time  base:  0.2  ms/div. 

Vertical  scale:  11.3  atm/div. 

Fig.  8  Pressure  records  showing  the  onset  of  single  head  spin  in 

ethylene-air  mixtures. 

a)  4.5%  C2H4  in  4.8  cm  diameter  tube.  Biomation  record: 
Time  base:  0.1  ms/div. 

Vertical  scale:  3.4  atm/div. 

b)  4.2%  in  4.8  cm  diameter  tube.  Biomation  record: 
Time  base:  0.1  ms/div. 

Vertical  scale:  6.8  atm/div. 

c)  4%  in  14.5  cm  diameter  tube.  Biomation  record: 

Time  base:  0.1  ms/div. 

Vertical  scale:  2.7  atm/div. 

d)  3.5%  in  14.5  cm  diameter  tube.  Biomation  record: 

Time  base:  0.1  ms/div. 

Vertical  scale:  6.8  atm/div. 


APPENDIX  IV 

On  the  Scaling  of  Blast  Waves  from  Fuel-Air  Explosives’ 


Paper  presented  at  the  Vlth  International  Symposium  on  Military 
Applications  of  Blast  Simulation,  Caiiors,  France,  June  1979. 
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In  a  conventional,  non-nuclear  weapon,  a  strong  blast  wave  of  high  energy 
density  is  generated  by  detonating  in  the  surrounding  atmosphere  a  concentrated 
charge  of  solid  (or  high)  explosive  (henceforth  H.E.).  Because  of  the  relative 
small  dimensions  of  the  concentrated  charge,  the  charge  dimensions,  its  geometry 
and  the  characteristic  time  for  the  energy  deposition  can  usually  be  ignored. 
Therefore,  blast  waves  generated  by  concentrated  charges  approach  ideal  point  blast 
behavior  at  relative  close  distances  from  the  source.  For  these  ideal  explosions, 
nearly  all  the  chemical  energy  released  by  the  detonation  process  can  be  considered 
to  be  transferred  to  the  shock  front  and  subsequently  dissipated  in  the  surround¬ 
ings  as  the  blast  wave  propagates.  Given  the  mass  of  the  solid  explosive  charge, 
the  blast  energy  Eq  is  thus  known  a  priori  for  ideal  explosions  and  a  characteris¬ 
tic  explosion  length  Rq  -  can  he  derived  to  scale  all  ideal  explosions. 

The  destructive  potential  of  conventional  weapons,  i.e.,  the  blast  damages,  can 
then  be  easily  estimated  from  a  standard  chart  giving  point  blast  overpressures 
Aps  and  static  and  dynamic  impulses  1^.  and  1^,  respectively,  with  the  scaled  dis¬ 
tance  R  /R  where  R  denotes  the  saock  radius^, 
so  s 

For  certain  military  applications,  blast  waves  generated  by  fuel-air  ex¬ 
plosives  (henceforth  FAC)  with  moderate  energy  densities  and  moderate  initial 
strengths  may  prove  more  effective.  In  an  FAE  weapon,  the  fuel  is  explosively  dis¬ 
seminated  in  the  atmosphere  to  form  an  explosive  vapor  cloud  and  then  detonated 
subsequently  by  one  or  more  initiating  ^barges.  Depending  on  the  rate  of  fuel  dis¬ 
persion  and  the  ignition  delay,  the  vapor  cloud  dimensions  (i.e.,  the  cloud  radius 
R  )  may  be  quite  large.  As  a  result  of  t ho  relative  large  dimension  of  the  cloud, 
vapor  cloud  explosions  generate  non-ideal  blast  waves  which  scale  according  to  the 

cloud  radius  R  in  the  near  field.  In  the  far  field  where  the  distances  are  large 
c 

compared  with  the  characteri s *  i c  dimetr.i  "i  of  l.h<>  cloud,  FAI  blast  waves  should 
dec, iv  accord  i  no  to  tho  Bothe-I  aiidati-l-'h  i  t  hum  asvmplotio  solution  as  in  the  case  of 
ideal  point  blast  waves,  hence  scale  according  to  the  explosion  length  R  . 


There- 
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fore,  FAC  weapons  require  two  characteristic  lengths  for  scaling  in  the  near  and 
far  fields,  respectively.  Whereas  the  cloud  radius  Rc  can  readily  he  estimated, 
the  explosion  length  in  FAE  weapons  can  no  longer  be  assessed  a  priori  like  in  a 
point  blast  explosion.  For  fuel-air  explosives,  the  actual  blast  energy  can  no 
longer  be  equated  to  the  total  combustion  energy  of  the  reactive  mixture  since  a 
significant  portion  remains  in  the  combustion  products  in  the  form  of  internal 
energy.  Furthermore,  because  of  the  large  dimensions  of  the  FAE  cloud,  the  blast 
decay  outside  the  cloud  can  be  drastically  affected  by  the  details  of  the  energy 
release  processes  inside  the  cloud.  In  other  words,  the  effective  blast  energy 
may  strongly  depend  on  the  various  possible  modes  of  combustion  of  the  cloud  (viz., 
detonation,  volumetric  explosion  or  deflagration).  The  effective  energy  should  be 
the  actual  work  done  by  the  expanding  interface  between  the  combustion  products 
and  the  surroundings.  The  blast  energy  cannot  be  assessed  a  priori  since  the  path 
of  the  interface  and  the  pressure  variation  on  the  interface  can  only  be  determined 
after  solving  for  the  blast  flow  structure  itself.  This  in  turn  requires  the  nume¬ 
rical  integration  of  the  gasdynamic  equations  with  the  i  ipropriate  initial  and 
boundary  conditions  of  the  problem.  Such  numerical  computations  ore  also  very  ex¬ 
pensive  because  good  accuracy  in  the  far  field  requires  the  use  of  very  fine  step 
sizes. 

To  estimate  the  blast  energies  from  fuel-air  and  oxygen  explosives,  the 

1  2 

Authors  have  used  the  simpler  Brinkley-Kirlcwood  theory  (henceforth  B-K)  and  match 

3 

the  far  field  experimental  overpressure  decay  with  distance  .  Good  agreement  was 
achieved  with  blast  energies  representing  approximately  20  and  25%  of  the  heat  re¬ 
leased  by  detonative  combustion  in  stoichiometric  hydrocarbon-air  and  oxygen  mix- 

lures,  respectively.  On  the  other  hand,  Fishlnirn  using  Oppenhein's  numerical  code 
r 

(viz.,  CLOUD  code)'  has  estimated  the  blast  enorqy  from  MAPP-air  detonation  at  ap¬ 
proximate!  v  31.?,"'  of  the  rombustion  energy.  The  CLOUD  code  does  not  calculate  the 
work  done  by  the  expanding  interface,  hut  the  blast,  energy  can  bo  estimated  from 
the  lirst  law  of  thermodynamics  and  the  interface  position  when  the  overpressure 
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at  ;  the  interface  returns  to  zero.  Since  MAPP-air  detonation  data^  are  very  simi¬ 
lar  to  the  hydrocarbon-air  data  used  in  Ref.  1,  it  is  difficult  to  explain  the 
discrepancy  in  both  estimates  for  fuel-air  mixtures,  the  CLOUD  code  blast  energy 
being  almost  twice  as  large  as  the  value  derived  from  the  B-K  theory. 

Furthermore,  blast  energy  estimates  from  volumetric  explosions  of  MAPP- 
air  and  hydrocarbon-air  mixtures  differ  only  by  a  factor  of  1.25.  For  MAPP-air, 
Fishburn's  blast  energy  estimate  of  33.6%  of  the  combustion  energy  is  almost  equal 
to  his  blast  energy  estimate  from  detonation  and  suggests  that  blast  energy  may  be 
independent  of  the  combustion  mode.  For  hydrocarbon-air  volumetric  explosion,  the 
Authors'  have  estimated  the  blast  energy  at  approximately  42%  of  the  combustion 
energy  by  assuming  an  isentropic  expansion  of  the  explosion  products  since  no  ex¬ 
perimental  data  were  available  to  fit  a  curve  of  overpressure  decay  with  distance 
as  performed  for  the  detonation  case.  This  value  which  constitutes  an  upper  bound 
for  the  work  done  by  the  expanding  interface  is  more  than  twice  the  Author's  blast 
energy  estimate  from  detonation  and  suggests  that  blast  energies  in  hydrocarbon-air 
mixtures  depend  on  the  mode  of  combustion  in  complete  disagreement  with  Fishburn's 
conclusion.  On  the  other  hand,  the  Authors'  blast  energy  estimate  from  volumetric 
explosion  in  hydrocarbon-oxygen  mixtures  of  2 ?.%  of  the  combustion  energy'  compares 
well  with  their  25%  estimate  from  detonation  and  tends  to  support  Fishburn's  con¬ 
clusion  regarding  the  independence  of  the  blast  energy  from  the  mode  of  combustion 
in  the  cloud. 

This  is  an  important  question  to  resolve  since  volumetric  explosion  may 
occur  if  partial  or  total  confinement  of  the  cloud  is  achieved  or  by  heating  uni¬ 
formly  the  cloud  to  its  auto-ignition  limit  by  radiation  or  shock  wave.  The  ini¬ 
tial  strength  of  the  blast  wave  from  volumetric  explosion  being  approximately  half 
II'1'  value  from  detonation  results  in  a  slower  decay  of  the  blast  strength,  hence, 
larger  dam, iocs  in  the  intermediate  field*' .  Therefore,  FAF  weapons  in  which  the  ex¬ 
plosive  vapor  cloud  is  vol  nine  t.ri  call''  oxolmied  may  be  quite  attractive  for  military 
appl icat ions. 
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,  Reliable  blast  energy  estimates  from  FAE  should  also  resolve  the  import¬ 

ant  problem  of  far  field  equivalency.  Far  field  equivalency  implies  that  ideal 
blast  wave s  from  H.E.  and  nonrdieal  blast  waves  from  FAE  are  equivalent  in  the 
far  field.  In  other  words,  on  the  basis  of  the  blast  energy  alone,  blast  overpres¬ 
sure  and  impulses  can  be  estimated  for  these  explosions  in  the  weak  shock  far  field 
regime. 

In  this  paper,  a  reliable  numerical  code  is  developed  to  compute  the 
blast  energy  from  vapor  cloud  explosions  for  far  field  scaling  of  non-ideal  blast 
waves  in  relation  with  FAE  weapons.  The  present  study  also  resolves  the  important 
questions  of  blast  energy  dependence  on  the  mode  of  combustion  in  the  cloud  and  far 
field  equivalency  of  H.E.  and  FAE  blast  waves  using  as  a  particular  example  the 
blast  waves  generated  by  centrally  ignited  detonations  and  volumetric  explosions  of 
spherical  clouds  of  stoichiometric  hydrocarbon-air  and  oxygen  mixtures. 


2 *  Vapor  Cloud  Explosion  Models 


2.i.  C-J  Detonation 


Assume  that  at  time  t  =  0,  a  C-J  detonation  is  initiated  at  the  center 
r  =  0  of  a  reactive  cloud  of  radius  Rc  and  propagates  at  a  constant  velocity  D  in 
the  reactive  mixture.  The  burnt  gases  behind  the  C-J  front  are  expanded  and  de¬ 
celerated  to  zero  velocity  by  a  centered  expansion  wave.  The  self-similar  flow 
field  behind  the  C-J  front  can  easily  be  computed  for  any  given  reactive  mixture 
with  prescribed  initial  conditions,  namely,  composition,  pressure,  temperature, 
specific  heat  ratios  of  unburnt  and  burnt  gases  and  combustion  energy^. 


When  the  detonation  front  reaches  the  cloud  edge  (r  =  Rc),  a  shock  wave 
is  transmitted  into  the  surrounding  atmosphere  and  an  exoansion  wave  propagates 
back  into  the  explosion  products.  The  subseouent  flow  structure  inside  and  outside 
the  cloud  is  no  longer  self-similar  and  can  only  be  described  by  solving  numeri¬ 
cally  the  time-dependent,  one-dimensional,  spherically  symmetric,  adiabatic 
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equations  for  compressible,  non-viscous  and  non-heatinq  conductinq,  perfect  qoses 
subjected  to  appropriate  boundary  conditions,  viz., 


dr¬ 

at 


+  dJ£ 

'dr 


o 


0) 


f>  =  fZT-  (4) 

where  o,  p,  T,  u  and  y  denote  the  density,  pressure,  temperature,  particle  velo¬ 
city  and  specific  heat  ratio,  respectively.  The  relevant  boundary  conditions  at 
the  shock  front  are  qiven  by  the  Rankine-Hugoniot  relations, 

--  (  Yo  +0/(To  -/  +  2 
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where  VJ  =  ,  Ms  donot.es  t ho  shock  f tach  number  and  subscripts  s  and  o  character 

ize  the  shock  front  and  the  initial  conditions  in  the  surrounding  atmosphere.  At 
the  center,  the  boundary  conditions  require  the  flow  to  be  at  rest,  i.e., 

u.(o,-ki  -  o  ,  l  m 


2. i i .  Volumetric  Explosion 

Assume  at  time  t  =  0,  a  volumetric  explosion  occurs  in  a  reactive  cloud 
of  radius  R  of  prescribed  composition  and  initial  conditions  as  in  2.i.  The  re- 

v« 

suiting  pressure  and  temperature  increases  can  easily  be  calculated  by  combining 
the  sceady-state,  conservation  of  mass,  momentum  and  energy  equations,  viz., 

A/4.  -  ~  CYb'iX'Idftt'O  (?) 

where 

^  -  /«  Q-lYtck.  (8) 

Subscripts  u  and  b  characterize  unburnt  and  burnt  gases,  respectively  and  Q  denotes 
the  combustion  energy  per  unit  mass  of  reactive  mixture.  In  deriving  Eg.  (7)  it 
has  been  assumed  that  the  molecular  weights  of  unburnt  and  burnt  gases  are  the 
same. 


The  subsequent  propagation  of  the  transmitted  shock  in  the  surrounding 
atmosphere  and  the  reflected  expansion  wave  in  the  explosion  products  are  obtained 
by  solving  numerically  the  unstead''  qasdynamic  equations  (Cos.  1-6)-. 

3.  Numerical  Analysis 


The  r. LOl 111  code  which  solves  t.ho  unsteady  qasdvnamc  equations  in 
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Lagrangian  coordinates  is  very  adequate  to  calculate  the  interface  trajectory  and 
properties  to  estimate  the  Mast  energy  from  vapor  cloud  explosions.  However,  the 
code  which  uses  the  von  Neimiann-Richtmeyer  artificial  viscosity  to  avoid  shock 
discontinuities  broadens  the  shock  front  and  also  creates  numerical  overshoots  and 
undershoots  which  must  be  removed  by  the  use  of  appropriate  smoothing  techniques. 

The  present  numerical  codes  uses  an  Fulerian,  finite-difference  algori¬ 
thm  which  combines  Van  Leer's  first-order  scheme7  with  Boris  et  al's  flux  correct- 

3-10 

ed  transport  (FCT)  technique  v  to  handle  accurately  steep  gradients  associated 
with  shock  propagation  and  prevents  large  numerical  overshoots  or  undershoots. 

In  Van  Leer's  first-order  scheme,  the  centered  difference  scheme  of  F.q. 
(1)  is  stabilized  by  a  diffusion  term,  viz., 

at 


+ 


(9) 


4  U) 


where  the  diffusion  coefficients  k  are  defined  as 


k?,± i&J  -  i  + 


and.* 


(10) 


-  86  - 


kma )  lu.a.ijrw^  on 

where  cn  and  un  are  the  local  sound  speed  and  particle  velocity  at  the  grid  point 
n,  respectively.  To  prevent  overshoots  near  sharp  gradients  the  time  step  At  is 
restricted  by  the  Courant-Friedrich-Lewy  stability  criterion 

-  6  llXm  j C rf\  -i  lUnO  (12) 

where  the  Courant  number  a  is  less  than  unity. 

The  FCT  technique  consists  of  three  finite-difference  operations:  a 
transport  and  a  diffusion  followed  by  an  antidiffusion.  In  the  FCT  technique, 
the  value  of  the  vector  F  in  Eq.  (1)  is  calculated  at  any  grid  point  n  at  time 
t  +  At  in  four  steps.  In  the  first  step  where  transport  and  diffusion  opera¬ 
tions  are  performed,  an  approximate  value  for  F^  at  time  t  +  At  is  calculated  by 
a  finite-difference  approximation  to  Eq.  (1)  at  time  t  +  At/2  and  an  additional 
diffusion  term,  viz., 

Fji+M)  =  F„(-k)  +  At  V&U  ■*- At/2.) 

d't 

~  ^  (13) 


The  diffusion  term  added  in  Eg.  (13)  is  an  error  with  respect  to  Eq.  (1)  and  is 
removed  in  the  fourth  step  by  adding  an  equal  but  opposite  antidiffusion 


wm*mea  ^rwdii 
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i,n  Eq.  (13). 


Some  provisional  fluxes  necessary  to  remove  the  three-point  diffusive 
error  term  are  calculated  in  the  second  step,  viz., 


04) 


These  fluxes  are  then  corrected  in  the  third  step  to  prevent  the  formation  of  new 
extrema  in  the  solution  or  the  amplification  of  ex i sting  extrema  when  used  to  per¬ 
form  in  the  fourth  step  the  antidiffusion  of  the  first  step.  The  antidiffusion 
fluxes  are  given  bv 

fo, 


faUi+W  *  sprt&n+l  -  m**l 


05) 


/777m  1 A 7i-£  J  I <$"n-bl  1 7  An-f-ljJ 


where 
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In  the  last  step,  antidiffusion  is  performed  to  obtain  the  value  of  the 
vector  F  at  time  t  +  At,  viz., 

M)  FJrU^l 

UC») 

+  4m- i 


This  technique  which  prevents  tho 


V 

-,N 

formation  of  extrema  associated  with 


the  diffusion  term  does  not  remove  valid  physical  extrema  arising  from  the  term 

f 

8F/3t  in  Co..  (1 ) . 

< 

Combining  Van  Leer's  scheme  with  the  FCT  technique,  i.e.,  substituting 
for  Eq.  (13),  the  first  step  in  FCT,  the  following  transport  and  diffusion  opera¬ 
tions, 

Et-t+M)  ^  F^tt)  +M  %6n(t+AF) 

m  dt 
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yields  less  expensive,  more  accurate  and  stable  solutions  to  some  yasdynamic  prob¬ 
lems  than  by  applying  each  method  separately.  Therefore,  the  Van  Leor-FGT  numeri¬ 
cal  method  has  been  used  to  compute  the  flow  structure  associated  with  vapor  cloud 
explosions . 

4 .  Re su Its  and  Discu ss_i on 

The  propagation  in  air  of  non-ideal  blast  waves  generated  by  centrally 
initiated  detonation  and/or  volumetric  explosion  of  spherical  clouds  of  stoichio¬ 
metric  propane-air  and  oxygen  mixtures  at  initial  pressure  and  temperature  of  1 
atmosphere  and  ?r>°C,  respective!'/,  have  boon  connoted  numerically  usinn  the 
Lulerian  code  described  in  Section  3.  Propane  was  chosen  as  the  fuel  because  it 
represent s  a  typical  hydrocarbon  and  furthermore  evoer iment al  blast  data  exist  for 
meaningful  comparison.  For  both  mixtures,  the  relevant  properties  corresponding 
to  both  combustion  modes  are  given  in  Table  1.  Because  of  the  differences 


in  the 


*i  * 


specific  heat  ratios  of  the  combustion  products  and  air,  the  specific  heat  ratio 
was  varied  smoothly  around  the  interface  from  the  explosion  product  value  to  the 

t 

air  value  to  prevent  the  formation  of  overshoots  around  the  interface.  The  inter¬ 
face  location  was  obtained  from  mass  conservation  of  the  explosion  products,  viz., 


R:W 


4xrpdr 
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where  denotes  the  instantaneous  location  of  the  interface  and  subscript  u  cha¬ 
racterizes  the  unburnt  mixture.  The  pressure  at  the  interface  was  calculated  by 
linear  interpolation.  Numerical  calculations  were  carried  out  to  blast  overpres¬ 
sures  of  approximately  0.1  atmosphere  and  the  FCT  diffusion  coefficient^  was  kept 
constant  at  a  value  of  1/8. 


Typical  pressure  profiles  for  a  propane-air  detonation  and  volumetric 
explosion, plotted  in  Figs.  1  and  2  with  respect  to  distance  scaled  according  to 
the  cloud  radius  Rj demonstrate  the  stability  of  the  numerical  technique  and  its 
ability  to  handle  sharp  discontinuities.  Similar  profiles  were  obtained  for  pro¬ 
pane-oxygen  mixtures  with  faster  decay  rates  associated  with  higher  initial  blast 
overpressures. 

Blast  energies  defined  as  the  work  done  by  the  expanding  interface  prior 
to  the  formation  of  the  negative  overpressure  phase  are  found  to  represent  apnroxi- 
mately  36,  34  and  27?;  of  the  combustion  energies  of  detonated  propane-air,  ethy¬ 
lene-air  and  propane-oxygen  mixtures,  respectively.  The  propane-oxygen  value  com¬ 
pares  well  with  the  Author's  previous  estimate  of  2D';  of  the  combustion  energy^. 

On  the  other  hand,  the  present  fuel-air  data  which  compares  well  with  Fishhurn's 
value  for  *1APP^  of  37.8";,  completely  disagrees  with  our  previous  estimate  of  al¬ 
most  20"' ^  ami  demonstrates  the  danger  of  using  fitting  techniques. 

For  volumetric  explosions  of  propane-air  and  oxvgen  mixtures,  the  blast 
energies  represent  approximately  34  and  278.  of  the?  respective  combustion  energies 
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of  these  mixtures  and  differ  from  the  corresponding  blast  encrqies  from  detonation 

i 

by  6  and  1",  respectively.  Therefore,  the  present  results  support  and  extend 
Fishburn's  conclusion  that  blast  energies  from  fuel-air  and  oxygen  mixtures  are 
almost  independent  of  the  mode  of  combustion  of  the  explosive  cloud  (i.e.,  detona¬ 
tion  or  volumetric  explosion).  The  larqer  percentage  of  combustion  energy  appear¬ 
ing  as  blast  energy  for  fuel-air  explosions  reflects  the  smallest  entropy  increase 
associated  with  fuel-air  combustion  (approximately  60"  less  than  for  fuel-oxygen 
mixtures).  Assuming  an  isentropic  expansion  of  the  combustion  products  to  ambient 
pressure,  the  percentages  of  combustion  energy  remaining  in  the  combustion  products 
as  internal  energy  have  been  estimated  at  about  57  and  70"  for  fuel -air  and  oxygen 
mixtures,  respectively,  in  good  agreement  with  the  corresponding  blast  energies. 

Non-ideal  blast  wave  overpressures  and  impulses  generated  by  fuel-air 
and  oxygen  explosions  plotted  in  Fig.  3  show  far  fi  Id  equivalency  when  scaled 
with  respect  to  the  explosion  length  R  .  The  B-K  theory  was  used  to  generate  the 
far  field  data.  Starting  at  5  cloud  radii,  the  B-K  data  wore  found  in  excellent 
agreement  with  the  numerical  data  which  were  computed  down  to  about  10  cloud  radii, 
except  for  the  propane-oxyqcn  static  impulses.  The  B-K  results  for  this  parameter 
are  20»  smaller  than  the  numerical  results  and  may  reflect  the  breakdown  of  the 
B-K  assumption  for  impulse  for  fuel -oxygen  mixtures. 

In  the  near-field,  blast  parameters  reflect  the  non-idealities  of  the 
explosion,  namely  the  mode  of  combustion  and  the  nature  of  the  oxidizer.  Volumet¬ 
ric  explosions  are  seen  to  be  more  destructive  in  the  intermediate  field  where  the 
blast  wave  decay  is  slower  than  the  decay  resulting  from  detonation  of  the  vapor 
cloud. 


Ideal  point  blast  parameters  also  plotted  in  Fig.  3  decay  faster  than 
the  non-ideal  data  in  t ho  intermediate  and  far  fields.  In  other  words,  st.rono 
blast  waves  of  high  energy  density  genera  tod  by  solid  explosives  are  less  destruc¬ 
tive  in  the  intermediate  and  tar  fields  than  the  moderate  strength  blast  waves 


associated  with  smaller  energy  density  explosions.  By  matching  the  far  field  beha¬ 
vior,  an  equivalent  blast  energy  ratio  E^  can  be  obtained.  In  the  present  study,  the 
the  overpressure  curves  have  been  matched  at  0.01  atmosphere  where  the  correspond¬ 
ing  scaled  distances  R/R  are  1.7  and  2.4  for  ideal  and  non-ideal  blast  waves,  res¬ 
pectively.  The  equivalent  blast  energy  ratio  E^  is  (1.7/2.4)3  =  2.0.  There¬ 
fore,  the  destructive  potential  of  a  non-ideal  blast  wave  from  fuel-air  or  oxygen 
detonation  or  volumetric  explosion  is  equivalent  to  that  of  an  ideal  point  blast 
explosion  with  nearly  three  times  the  explosion  enerqy.  This  is  a  consequence  of 
the  larger  dissipation  of  energy  in  point  blast  with  high  initial  overpressure.  Ry 
the  time  the  point  blast  wave  has  decayed  to  an  acoustic  wave,  almost  no  energy  is 
left  in  the  N-wave.  In  contrast,  non-ideal  blast  waves  of  moderate  initial  over¬ 
pressure  dissipate  less  energy  as  they  propagate,  hence  still  retain  in  the  far- 
field  a  significant  portion  of  their  energy. 

Using  the  effective  blast  wave  energy  for  scaling  non-ideal  blast  waves, 
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theoretical  blast  parameters  are  compared  in  Fig.  4  with  experimental  data  .  The 
latter  were  obtained  in  spherical  balloons  of  diameters  ranging  from  0.7  to  6 
meters,  centrally  ignited.  The  data  are  representative  of  detonations  in  several 
hydrocarbon-air  and  oxygon  mixtures  and  show  no  dependence  on  the  fuel  and  a  slight 
dependence  on  the  nature  of  the  oxidizer  (air  or  pure  oxygen).  Empirical  formulae 
derived  in  Ref.  3  to  fit  t.he  experimental  data  have  also  been  plotted  in  Fig.  4. 

The  measurements  were  carried  out  to  distances  of  50  cloud  radii.  Theoretical 
overpressures  in  the  far  field  differ  from  the  experimental  data  at  most  by  10%, 
well  within  the  5-10."  accuracy  of  the  measurements.  The  corresponding  static  im¬ 
pulses  show  less  agreement  (20%  greater  than  the  experimental  data)  in  view  of  the 
difficulties  associated  with  static  impulse  measurements. 

5.  Conclusions 

A  Lulerian  numerical  code  usinq  a  f lux-corrected  transport  technique 
capable  of  handling  accurately  sharp  discontinuities  1 ike  shock  waves  has  been 


developed  to  c  .nputo  the  effective  energy  of  non-ideal  blast  waves  generated  by 
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FAF.  weapons.  For  most  hydrocarbon-air  mixtures,  it  is  found  that  almost  271  of 
the  combustion  energy  actually  goes  to  the  blast  wave  and  this  fraction  decreases 
for  mixtures  with  higher  energy  densities  (eg.,  about  27"  for  fuel-oxygen  mixtures). 
Furthermore,  the  effective  blast  energy  of  fuel-air  and  oxygen  mixtures  exhibit  no 
significant  dependence  on  the  mode  of  combustion  of  the  reactive  mixture  whether 
detonation  or  volumetric  explosion.  On  the  basis  of  the  effective  blast  energy, 
non-ideal  blast  waves  from  fuel-air  and  oxygen  explosives  are  all  equivalent  in 
the  far  field  and  the  blast  parameters  (eg.,  overpressure,  static  and  dynamic  im¬ 
pulses)  scale  according  to  the  explosion  length.  However,  far  field  equivalence 
for  ideal  and  non-ideal  blast  waves  does  not  exist.  In  other  words,  far  field 
blast  parameters  from  ideal  and  non-ideal  explosions  arc  not  equal  at  the  same 
scaled  distance.  By  matching  the  blast  overpressure  curves  in  the  far  field,  the 
equivalent  point  source  energy  for  the  non-ideal  blast  wave  from  a  fuel-air  or 
oxygen  explosive  weapon  is  three  times  greater  than  the  energy  from  a  point  blast. 
As  a  result,  the  destructive  potential  of  non-ideal  blast  waves  from  FAF  weapons 
is  found  to  be  superior  to  that  of  conventional  II.F.  weapons  in  the  far  field. 
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Table  1 


Volumetric 
Detonation  Explosion 


Mi xture  ] 

Fuel,  % 

Yu 

% 

( a  tm ) 

Yb 

APV 
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Yb 

Q 

(kcal/kg) 

CaHg-Air  i 

4. 

1.37 

17.27 

1.28 

8.34  ! 

1.26 

668 

C2H4-A1r 

6.53 

1.38 

17.74 

1.27 

8.32 

1.25 

720 

C3H8_02 

16.67 

j  1.29 

35.2 

1.13 

17.2 

1.136 

2400 

Pressure  profiles  from  stoichiometric  propane-air  volumetric  explosion 
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